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FEDERAL  AVIATION  AIMINISTRATION 
SYSTEMS  RESEARCH  AND  DEVELOFMEWT  SERVICE 
SPECTRUM  MANAGEMENT  STAFF 


STATEMENT  OF  MISSION 


uiiiiSion  of  the  Spectrum  Management  Staff  is  to  assist  the  Department  of 
State,  Office  of  Telecommunications  Policy,  and  the  Federal  Communications 
Corrnii.SfciJ.on  in  assuring  the  FAA's  and  the  nation's  aviation  interests  with 
sufficient  protected  electromagnetic  telecommunications  resources  throughout 
the  world  to  provide  for  the  safe  conduct  of  aeronautical  flight  by  fostering 
ef '’oc t.’ ve  and  efficient  use  of  a natural  resource  - the  electromagnetic  radio 
frequency  spectrum. 

This  objective  is  achieved  through  the  following  services: 

^ Planning  and  defending  the  acquisition  and  retention  of 
sufficient  radio  frequency  spectrum  to  support  the  aero- 
nautical interests  of  the  nation,  at  home  and  abroad,  and 
spectrum  standardization  for  the  world's  aviation  community. 

• Providing  research,  analysis,  engineering,  and  evaluation 
in  the  development  of  spectrum  related  policy,  planning, 
standards,  criteria,  measurement  equipment,  and  measure- 
ment techniques. 


• Conducting  electromagnetic  compatibility  analyses  to 
determine  intra/lnter-system  viability  and  design 
parameters,  to  assure  certification  of  adequate  spectrum 
to  support  system  operational  use  and  projected  growth 
patterns,  to  defend  aeronautical  services  spectrum  from 
encroachment  by  others,  and  to  provide  for  the  efficient 
use  of  the  aeronautical  spectrum. 

• Developing  automated  frequency  selection,  computer  programs/ 
routines  to  provide  frequency  planning,  frequency  assign- 
ment, and  spectrum  analysis  capabilities  in  the  spectrum 
supporting  the  National  Airspace  System. 


Providing  spectrum  management  consultation,  assistance, 
and  guidance  to  all  aviation  interests,  users,  and 


international . 
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OBJECTIVE  MEASUREMENT  OF  VOICE 
CHANNEL  INTELLIGIBILITY 

K.J.  Gamauf 
W.J.  Hartman* 


Following  the  results  of  a feasibility 
study  (Hartman  and  Boll,  1976)  an  objective 
intelligibility  measure  is  developed  using 
a large  data  base  consisting  of  8-50  word 
phonetically  balanced  word  groups  with  twelve 
different  kinds  of  distortion.  Justifi- 
cation for  the  use  of  this  particular  measure 
is  included,  with  mathematical  derivations 
and  physical  interpretations. 

A discussion  of  the  feasibility  of  a 
hardware  implementation  of  the  software 
developed  here  is  also  included. 

Key  words:  intelligibility  measurements; 

linear  predictive  coding;  voice 
systems . 


1.  INTRODUCTION 

There  has  long  been  a need  for  an  inexpensive,  reliable,  and 
efficient  method  to  evaluate  the  quality  of  speech  sent  over  voice 
communication  channels.  Few  voice  communication  systems  today 
are  judged  by  the  quality  and  intelligibility  of  the  speech 
received  by  the  listener.  Instead,  system  performance  is 
generally  specified  by  some  engineering  parauneter,  such  as  the 
signal-to-noise  ratio  of  the  receiver  output. 

The  most  common  procedure  for  determining  the  intelligibility 
of  a voice  channel  is  a subjective  method  that  involves  trained 
speakers  and  listener  panels  that  directly  score  the  percentage 
of  speech  that  is  intelligible.  These  schemes  have  the  desirable 
property  that  they  produce  repeatable  results.  Unfortunately, 
subjective  scoring  methods  are  expensive  and  time  consuming 
and  as  a result,  are  not  widely  used.  What  is  needed  is  an 

*The  authors  are  with  the  Institute  for  Telecommunication  Sciences, 
Office  of  Telecommunications,  U.S.  Department  of  Commerce, 

Boulder,  CO  80302. 


inexpensive  and  efficiently  applied  objective  evaluation  of 
speech  intelligibility  that  is  comparable  to  subjective  methods. 

This  paper  develops  a method  of  obtaining  an  objective 
intelligibility  measure  that  gives  good  results  for  speech  sent 
through  both  analog  and  digital,  noise-corrupted  communication 
channels.  The  distortion  measure  is  obtained  using  Linear  Pre- 
dictive Coding  (LPC) , a mathematical  technique  widely  known  for 
its  application  to  the  analysis  and  synthesis  of  speech.  The 
feasibility  of  using  LPC  to  develop  an  objective  intelligibility 
measure  has  been  demonstrated  by  Hartman  and  Boll  (1976) . 

A.  The  Articulation  Index 

A well  known  objective  measure  that  is  used  for  voice  communi- 
cation channel  evaluation  is  the  Articulation  Index  (AI)  (ANSI, 
1969;  Kryter,  1962) . The  AI  is  a physical  measurement  that  is 
highly  correlated  to  speech  intelligibility  under  certain 
conditions.  The  AI  is  obtained  by  evaluating  the  signal-to-noise 
power  ratios  in  20  specific  frequency  bands.  These  power  ratios 
are  then  summed  and  normalized  to  give  a score  between  zero  and 
one.  An  AI  equal  to  one  signifies  perfect  intelligibility,  while 
a value  of  zero  represents  a complete  lack  of  intelligibility. 

The  AI  can  be  computed  using  a general  purpose  computer  once  the 
required  spectral  data  are  obtained  from  the  speech. 

An  automated  technique  to  obtain  the  AI  is  achieved  through 
the  Speech  Communication  Index  Meter  (SCIM) (Kryter  and  Ball, 

1964) , which  uses  9 frequency  bands,  instead  of  20,  to  obtain  a 
modified  AI . The  SCIM  system  can  be  used  to  perform  on-line 
measurements  and  has  been  found  to  be  highly  correlated  with  the 
standard  20  band  AI.  The  SCIM  system's  AI  can  also  be  directly 
related  to  speech  intelligibility  as  long  as  the  noise  present  is 
generally  additive  white  gaussion  noise.  Clipping  of  the  speech 
by  the  voice  channel  or  noise  that  is  intermittent  or  colored 
distorts  the  AI,  and  a correction  factor  must  be  employed. 
Multiplicative  noise  requires  a complete  recalibration  of  the 
SCIM  system.  Therefore,  the  type  of  noise  or  distortion  present 
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in  the  voice  channel  must  be  knovm  in  order  to  obtain  accurate 
results.  When  digital  voice  systems  were  tested  on  the  SCIN 
scheme,  very  poor  estimates  of  speech  intelligibility  were 
obtained.  This  was  true  even  when  quantization  noise  was  th6 
only  distortion  present.  Reliable  correction  factors  for  digital 
voice  systems  have  as  yet  not  been  found  to  compensate  for  this 
poor  performance.  A more  detailed  description  of  the  AI  and  the 
SCIM  system  of  objective  voice  channel  evaluation  can  be  found  in 
the  work  of  Hubbard  and  Hartman  (1974) . 

B.  Chapter  Summaries 

Chapter  2 gives  a brief  discussion  of  linear  prediction  of 
the  speech  waveform.  Chapter  3 describes  the  analog  voice  tapes 
used  in  this  study.  Chapter  4 describes  the  data  processing  of 
the  analog  voice  tapes  in  order  to  obtain  the  LPC  information 
from  the  speech.  Chapter  5 discusses  the  distance  or  distortion 
measures  that  were  used  to  predict  objectively  the  intelligibility 
of  noise  corrupted  speech.  A comparison  between  objective 
distance  measures  and  subjective  intelligibility  scores  is  given. 
Chapter  6 gives  a block  diagram  for  hardware  implementation  of 
the  scoring  techniques.  Appendix  A contains  listings  of  the 
computer  programs  used  to  obtain  the  numerical  results  found  in 
this  report. 

2.  LINEAR  PREDICTIVE  CODING 

Linear  predictive  coding  (LPC)  has  long  been  used  in 
communication  theory.  More  recently,  it  has  found  applications 
in  speech  analysis  and  synthesis,,  speaker  identification,  and 
word  recognition,  to  name  just  a few  new  areas.  In  this  study, 

LPC  is  used  to  develop  an  objective  intelligibility  measure  of 
speech  corrupted  by  noise. 


A.  Linear  Prediction 

LPC  models  the  vocal  tract  as  an  all-pole  digital  filter  and 
estimates  the  filter  parameters  (predictor  coefficients)  using 
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the  time  domain  speech  waveform  itself,  rather  than  the  waveform's 
short-term  frequency  spectrum.  This  makes  LPC  a relatively 
efficient  method  for  encoding  speech  compared  to  frequency  domain 
techniques. 

The  vocal  tract  is  assumed  to  be  modeled  as  a discrete, 
time-varying  filter  with  parameters  changing  slowly  enough  so  that 
they  can  be  considered  fixed  over  a specified  time  interval. 

Hence,  the  vocal  tract  can  be  approximated  by  a series  of 
stationary  shapes.  Atal  and  Hanauer  (1971)  have  shown  that  this 
all  pole  model  can  account  for  the  glottal  volume  flow  and 
radiation  of  sound  from  the  mouth  in  addition  to  vocal  tract 
sounds . 

The  transfer  function,  H(z),  used  to  describe  the  digital 
model  over  each  analysis  frame  is  given  by 


H(z)  = 


1-  E 

i=l 


-i 


(2.1) 


for  a model  with  P poles. 

The  time  sequence  corresponding  to  the  output  of  the 
recursive  filter  can  be  written  as 

P 


S = 
n 


^i  ^n-i 


+ 6. 


n 


n = 0,  1, 


(2.2) 


i=l 


where  (a^^)  are  the  predictor  coefficients  that  completely 
describe  the  characteristics  of  the  filter  and  (6^^)  is  the  driving 
function  or  input  to  the  filter. 

While  there  have  been  several  formulations  for  the  estimation 
of  the  linear  prediction  coefficients,  two  least  squares  methods 
have  become  prominent,  the  autocorrelation  method  and  the  co- 
variance  method.  The  autocorrelation  method,  which  will  be 
justified  in  Chapter  4,  was  chosen  for  use  in  this  study  and 
further  discussion  will  be  restricted  to  that  scheme.  The  auto- 
correlation method  can  be  considered  as  estimating  the  filter 
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coefficients  by  approximating  the  spectrum  of  the  speech 
waveform  by  an  all-pole  model. 

The  portion  of  the  signal  to  be  analyzed  is  first  multiplied 
by  a finite  window  of  length  N,  changing  the  signal  to 


Sn  = 


[Windowed  speech  samples,  o£n£N-i 


o,  n<o  and  n>N 


(2.3) 


Using  this  windowed  signal,  the  prediction  error  sequence  is 
defined  as 


'r.  = a.  S . 

n n 1 n-i 


n = 0,  1,  . . 


(2.4) 


and  the  total  squared  error  is  then 


N-l+P 


E e ^ = S ^ + 

n o 


E - E ’i 


(2.5) 


The  predictor  coefficients  are  selected  so  as  to  minimize 
the  total  squared  error.  This  is  accomplished  by  setting  the 
partial  derivative  of  the  total  squared  error  with  respect  to 
each  predictor  coefficient  equal  to  zero.  The  system  of  equations 
that  results  is 


E^u- 


i-k|  = ■'i  i - 1.  2 P 


(2.6) 


where 


N-l-li| 

i = R"  ^n  ^n+|i| 


(=  0,  1,  2,  . . . , P (2.7) 


are  the  normalized  short-term  autocorrelation  values  of  the 

speech  signal  and 
N-1 


R = ^ S “ 

o n 


(2.8) 


1 


is  the  normalization  factor  for  these  values.  The  normalized 
total  squared  error  can  be  defined  by  making  use  of  equations 
(2.5)  and  (2.6),  yielding 


E . = 1 

min 


- ■ 


(2.9) 


i=l 


The  predictor  coefficients  are  obtained  by  inverting  a 
positive  definite  Toeplitz  matrix 


i-k|l  . V _ 1 9 

'Ji,  k-1,  2,  . . •/ 


(2.10) 


This  system  of  equations  can  be  solved  by  using  Levinson's 
recursion  method,  which  will  be  expounded  upon  further  in 
Chapter  4.  The  Toeplitz  matrix  is  sometimes  called  the  auto- 
correlation matrix,  and  the  coefficients  obtained  from  this 
linear  system  result  in  a recursive  filter,  H(z),  which  is 
guaranteed  to  be  stable,  (all  of  its  poles  lie  inside  the  unit 
circle) , as  shown  by  Grenander  and  Szego  (1958) . 

B.  Spectral  Approximation 

Further  insight  can  be  gained  by  looking  at  the  frequency 
domain  approximation  to  the  above  system.  Taking  the  z-transform 
of  equation  (2.4),  one  obtains 


E(z)  = S(z)  [H(z)  ] 


-1 


(2.11) 


where  H(z)  is  defined  in  equation  (2.1)  and  E(z)  and  S(z)  are  the 
z-transforms  of  E^  and  respectively.  Rearranging,  equation 
(2.11)  can  be  written  as 


S(z)  = E(z)  H (z)  . 


(2.12) 


Minimizing  the  total  squared  prediction  error  is  equivalent 
to  approximating  the  error  sequence,  (e^) , by 


I 


®n 


lA  , n = o 


O , n o 


(2.13) 


in  a least  squares  sense.  This  implies  that  E(z)  is  being 
approximated  by  the  function  A,  a constant,  and  S(z}  is  being 
approximated  by  a spectrvim  corresponding  to  an  all-pole  transfer 
function,  i.e.. 


E(z)  = A 


(2.14) 


S{z)  = E(z)  H(z)  


-Z 


(2.15) 


The  value  of  A is  determined  by  the  application  of  energy  con- 
servation between  and  e^.  Using  equations  (2.9)  and  (2.13), 
one  obtains 


A^  = E . = 1 

min 


- Z 


(2.16) 


thereby  showing  that  A^  is  equal  to  the  minimum  total  squared 
error  of  the  system. 

This  approach  of  estimating  filter  coefficients  so  as  to 
minimize  the  energy  of  the  output  of  the  inverse  of  a system 
driver  by  its  impulse  response  is  sometimes  called  deconvolution 
or  inverse  filtering.  Considerable  work  has  been  done  in  this 
area  of  linear  prediction  of  speech  in  the  past  few  years  and 
many  good  references  are  available  that  give  more  detailed 
discussions  of  this  subject.  Some  particularly  good  ones  are 


Markel  and  Gray  (1976  and  1973);  Makhoul  (1975  and  1973); 
Makhoul  and  Wolf  .(1972);  and  Boll  (1973). 


3.  THE  VOICE  TAPES 


il 

if 

i In  order  to  develop  an  objective  intelligibility  measure  for 

corrupted  speech,  a comparison  must  be  performed  between  the 
j distorted  speech  and  the  original  noise  free  speech.  A subjective 

intelligibility  measure  of  the  distorted  speech  must  also  be 
available  in  order  to  judge  the  quality  of  the  objective  measure 
||  being  used.  Both  of  these  requirements  are  met  by  first  making  a 

noise  free  master  tape  of  pre-selected  speech  and  then  sending  it 
; through  voice  communication  channels  to  be  tested  and  making  a 

■ recording  of  the  speech  at  the  channel  output.  This  recording 

I can  be  subjectively  scored  for  intelligibility  and  also  compared 

f to  the  original  speech  by  some  mathematical  technique  to  obtain 

5 an  objective  measure. 


I 

'f 

\ 

I 


I 
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A.  Description  of  Voice  Tapes 

The  pre-selected  speech  to  be  sent  over  a voice  channel 
for  intelligibility  scoring  are  phonetically  balanced  (PB)  groups 
of  isolated  words  as  opposed  to  complete  sentences  or  nonsense 

4 

syllables.  These/PB  words  were  used  because  subjective  scores 
have  been  shown  'to  be  repeatable,  which  is  a necessary  criterion 
for  this  study  because  the  objective  measure  will  be  repeatable. 
Eight  PB  word  groups,  each  containing  fifty  isolated  words  were 
selected  as  the  test  speech.  A list  of  the  fifty  words  in  each 
word  group  is  given  in  Tables  3-1  through  3-8,  with  their 
designated  word  group  numbers. 

An  analog  tape  containing  all  eight  word  groups  and  using 
both  male  and  female  trained  speakers  was  obtained  from  the  Army 
Electronic  Proving  Ground  Electromagnetic  Environment  Test 
Facility  at  Fort  Huachuca,  Arizona.  From  this  tape,  a master 
analog  tape  was  made  that  would  be  sent  over  voice  channels  and 
later  compared  with  the  recorded  output  of  the  channel.  In  order 
to  perform  this  comparison,  the  two  tapes  would  have  to  be  aligned, 
which  meant  synchronization  information  must  be  included  on  the 
master  tape  before  being  sent  across  the  voice  channel.  Because 
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Table  3-1 
PB  Word  Group  361 


STAB 

TUCK 

DRAPE 

PITCH 

INK 

AID 

KIND 

STRESS 

TURN 

DROOP 

PUMP 

SUIT 

BARGE 

KNEE 

DUB 

WIELD 

ROCK 

BOOK 

THOU 

LAY 

FIFTH 

ROGUE 

CHEESE 

LEASH 

FRIGHT 


26. 

27. 

28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48. 

49. 

50. 


RUG 

CLIFF 

LOUSE 

G7^ 

RYE 

SANG 

CLOSED 

THREE 

MAP 

GAS 

SHEEP 

CREWS 

THRESH 

NAP 

HAD 

SHEIK 

TIRE 

DAME 

NEXT 

HASH 

SOAR 

TON 

DIN 

PART 

HOSE 


I 


Table  3-2 


PB  Word  Group  312 

1. 

JAB 

26. 

DIP 

2. 

ARC 

27. 

URGE 

3. 

JAUNT 

28. 

MOUTH 

4. 

ARM 

29. 

NET 

5. 

SHOP 

30. 

WAVE 

6. 

BEAM 

31. 

FINE 

7. 

KIT 

32. 

PURSE 

8. 

BLISS 

33. 

GOAT 

9. 

SPRIG 

34. 

HOG 

10. 

LAG 

35. 

RISK 

11. 

CHUNK 

36. 

DOUBT 

12. 

LATCH 

37. 

PUNK 

13. 

CODE 

38. 

DRAKE 

14. 

LOW 

39. 

WOOD 

15. 

TAB 

40. 

FEEL 

16. 

SHOT 

41. 

PROD 

17. 

SIGN 

42. 

FRISK 

18. 

CRUTCH 

43. 

DULL 

19. 

SAP 

44. 

MOST 

20. 

LOSS 

45. 

FUDGE 

21. 

CLASH 

46. 

POND 

22. 

SNOW 

47. 

HAVE 

23. 

CRY 

48. 

REEF 

24. 

SPY 

49. 

PROBE 

• 

in 

STIFF 

50. 

RICE 

10 


Table  3-3 
PB  Word  Group  291 


26 . NUTS 


1.  ARCH 


27.  WIPE 


29.  WITH 


5.  BIRTH 


31.  NERVE 


6 . SMART 

7.  SPUD 

8 . CLUB 


FLUFF 


NOOSE 


35.  FUME 

36.  WEAK 

37.  FUSE 


10 . CROWD 

11.  THAN 


WILD 


13.  THANK 


40.  PHONE 


15.  THRONE 

16 . CARVE 


42.  HOOF 


17.  TOAD 


19 . CHESS 

20 . TROOP 


21.  CHEST 

22.  BOOST 

23.  CLOWN 


46 . ITCH 

47.  ROOT 

48.  GRACE 


24.  DITCH 


50 . RUDE 


Table  3-4 
PB  Word  Group  265 


1. 

AS 

2. 

BEST 

3. 

EAT 

4. 

THUS 

5. 

EYES 

6. 

SCAN 

7. 

COB 

8. 

FALL 

9. 

DAD 

10. 

ODE 

11. 

SHANK 

12. 

MASH 

13. 

HITCH 

14. 

ROUGH 

15. 

FEE 

16. 

CHART 

17. 

WASP 

18. 

HULL 

19. 

TONGUE 

• 

o 

CM 

PUN 

21. 

REAP 

22. 

PUS 

23. 

BADGE 

24. 

DEEP 

25. 

SLOUCH 

26. 

CLOTH 

27. 

GROPE 

28. 

KEPT 

29. 

RAY 

30. 

FORGE 

31. 

CLOTHES 

32. 

ROOMS 

33. 

LAG 

34. 

THIGH 

35. 

WAIT 

36. 

WIFE 

37. 

JAG 

38. 

NIGH 

39. 

CRIB 

40. 

PRIG 

41. 

FLOP 

42. 

SUP 

43. 

GAGE 

44. 

WRIT 

45. 

PRIME 

46. 

FOWL 

47. 

BOG 

48. 

GAP 

49. 

FLICK 

50. 

RAISE 

12 
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Table  3-5 
PB  Word  Group  275 


1. 

AM 

26. 

SLEDGE 

2. 

GRADE 

27. 

RANGE 

3. 

GASP 

28. 

WOO 

4. 

MOTE 

29. 

DOPE 

5. 

MUD 

30. 

FLING 

6. 

BY 

31. 

NINE 

7. 

PHASE 

32. 

SCOUT 

8. 

RASH 

33. 

OFF 

9. 

RICH 

34. 

PIG 

10. 

POUNCE 

35. 

FORT 

11. 

SHAFT 

36. 

WOE 

12. 

ROAR 

37. 

CHOP 

13. 

ACT 

38. 

PLOD 

14. 

AIM 

39. 

KNIT 

15. 

HIM 

40. 

WHIFF 

16. 

COAST  • 

41. 

PENT 

17. 

DOSE 

42. 

THOUGH 

18. 

BUT 

43. 

JUG 

19. 

SOUTH 

44. 

SNIFF 

20. 

SIEGE 

45. 

QUIZ 

21. 

DWARF 

46. 

GUN 

22. 

FAKE 

47. 

COOK 

23. 

CUT 

48. 

SAG 

24. 

COMES 

49. 

WIRE 

25. 

SIN 

50. 

RAID 

13 
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Table  3-6 
PB  Word  Group  305 


1. 

STAFF 

26. 

TREE 

2. 

BASH 

27. 

GOOSE 

3. 

HAT 

28. 

PAGE 

4. 

WADE 

29. 

MAZE 

5. 

CHAMP 

30. 

FLIGHT 

6. 

ETCH 

31. 

PINK 

7. 

SLUG 

32. 

BUG 

8. 

CHANCE 

33. 

RAPE 

9. 

WAKE 

34. 

EARS 

10. 

VALVE 

35. 

SCRUB 

11. 

YOUTH 

36. 

COW 

12. 

FLAUNT 

37. 

TAG 

13. 

RUSH 

38. 

JAY 

14. 

GULL 

39. 

VOID 

15. 

DAUB 

40. 

EARTH 

16. 

REAL 

41. 

THOSE 

17. 

AIL 

42. 

LAP 

18. 

PUT 

43. 

SNIPE 

19. 

NUDGE 

44. 

FIR 

20. 

BACK 

45. 

CLOTHE 

21. 

PLUS 

46. 

MOPE 

22. 

BOB 

47. 

CORD 

23. 

THUG 

48. 

RIP 

• 

CM 

CUE 

49. 

HURT 

25. 

LINE 

50. 

FORCE 

1 
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Table  3-7 


PB  Word  Group  214 


1. 

TOE 

26. 

HID 

2. 

ARE 

27. 

SUCH 

3. 

RUB 

28. 

CRASH 

4. 

GROVE 

/ 29. 

BOX 

5. 

PANTS 

30. 

THERE 

6. 

DEATH 

31. 

END 

7. 

BAD 

32. 

MANGE 

8. 

PAN 

33. 

PLUSH 

9. 

USE 

34. 

IS 

10. 

SLIP 

35. 

FORD 

11. 

BASK 

36. 

HUNT 

12. 

FRAUD 

37. 

RAG 

13. 

NOT 

38. 

FEAST 

14. 

DEED 

39. 

NO 

15. 

SMILE 

40. 

CLOVE 

16. 

DISH 

41. 

FERN 

17. 

RISE 

42. 

PILE 

18. 

FUSS 

43. 

STRIFE 

19. 

WHEAT 

44. 

CANE 

20. 

DIKE 

45. 

FOLK 

21. 

PEST 

46. 

RAT 

22. 

CREED 

47. 

CLEANSE 

23. 

HEAP 

48. 

THEN 

24. 

BAR 

49. 

RIDE 

25. 

NOOK 

50. 

HIVE 
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Table  3-8 
PB  Word  Group  283 


1. 

US 

26. 

BIND 

2. 

SHACK 

27. 

CHEW 

3. 

CRACK 

28. 

WHEEZE 

4. 

CHANT 

29. 

FREAK 

5. 

YEAST 

30. 

PINT 

6. 

ASK 

31. 

GUESS 

7. 

EASE 

32. 

QUEEN 

8. 

REST 

33. 

CLOD 

9. 

JELL 

34. 

LOOK 

10. 

BOLT 

35. 

FRONT 

11. 

KILL 

36. 

NIGHT 

• 

CM 

LICK 

37. 

WIG 

13. 

CALF 

38. 

ROPE 

14. 

CATCH 

39. 

DAY 

15. 

TILL 

40. 

RHYME 

16. 

EACH 

41. 

SLIDE 

17. 

ROT 

42. 

FROCK 

18. 

ROLL 

43. 

LEFT 

19. 

BID 

44. 

FOOD 

20. 

COD 

45. 

SPICE 

21. 

DEUCE 

46. 

BORED 

22. 

DUMB 

47. 

THIS 

23. 

FAD 

48. 

THREAD 

24. 

HUM 

49. 

FORTH 

25. 

ROD 

50. 

FLIP 
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the  tapes  would  be  processed  in  a digital  state,  the  alignment 
procedure  would  also  have  to  work  in  a digital  format.  It  was 
found  that  a shift  of  plus  or  minus  10  samples  of  a 256  sample 
analysis  window  caused  the  predictor  coefficients  to  vary  less 
than  0.1%  in  all  cases.  Therefore,  the  synchronization  procedure 
to  be  used  was  required  to  align  two  segments  of  digitized 
speech  to  within  10  samples. 

It  should  be  noted  that  the  bound  on  the  variation  of  the 
predictor  coefficients  cannot  be  translated  into  a bound  on  the 
distortion  measures  described  in  Chapter  5.  However,  the 
alignment  method  described  in  the  next  section  was  tested 
extensively,  and  never  produced  a variation  in  the  distortion 
measures  larger  than  that  produced  by  the  normal  round  off  error. 

B.  Voice  Tape  Alignment 

A synchronization  procedure  that  was  found  to  meet  the 
required  10  sample  variation  specification,  made  use  of  a binary 
pseudo  noise  (PN)  sequence.  The  binary  PN  sequence  was  sent 
through  a phase-continuous  frequency  shift  keying  modem  using  the 
two  frequencies  1.2  kHz  and  2.2  kHz.  Several  different  length 
binary  PN  signals  and  modem  bit  rates  were  tested  to  determine 
the  best  combination  for  alignment  capability.  The  test  consisted 
of  cross-correlating  the  PN  sequences  under  different  noise  and 
distortion  conditions  and  looking  for  an  impulse  like  correlation 
function.  A further  requirement  was  to  have  the  PN  sequence  as 
short  as  possible.  It  was  found  that  a length  127  binary  PN 
sequence  sent  through  the  FSK  modem  operating  at  a 635  Hz 
bit  rate  followed  by  a low  pass  filter  with  a cutoff  of  2.5  kHz 
met  all  the  requirements  necessary  to  insure  the  alignment  of 
two  PN  sequences  distorted  by  noise.  The  low  pass  filter  was  used 
to  make  certain  that  the  frequency  spectrum  of  the  PN  signal  was 
in  the  range  required  for  input  to  most  voice  systems.  A PN 
signal  was  then  placed  before  each  word  and  after  the  last  word 
of  all  eight  word  groups  thereby  creating  the  master  analog  tape 
with  alignment  capabilities. 
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In  order  to  align  a distorted  tape  with  the  master  tape,  the 
location  of  all  the  PN  sequences  and  words  on  the  digitized 
master  tape  had  to  be  known.  This  was  done  by  blocking  the 
quantized  samples  into  125  sample  records  and  computing  the  mean 
and  standard  deviation  (SD)  for  each  record.  The  SD  was  used  as 
an  energy  criterion  to  determine  the  midpoints  of  the  PN  sequences 
and  words  and  the  length  of  each  word.  The  distances  between  the 
midpoints  of  each  PN  sequence  and  the  word  following  it  were  then 
determined.  The  corresponding  midpoints  of  the  words  from  a 
distorted  tape  are  now  all  that  remains  to  be  found. 

Each  of  the  eight  word  groups  on  a digitized  tape  made  up 
one  file  and  corresponding  files  between  the  master  and  distorted 
tape  were  aligned  independently  from  the  other  seven  sets  of 
files.  Using  the  SD  energy  criterion,  the  midpoints  of  the  first 
and  last  PN  sequences  of  the  distorted  word  groups  are  estimated. 
The  cross-correlation  between  these  PN  sequences  and  the  corres-  ^ 
ponding  ones  from  the  master  tape  are  then  computed,  thereby 
obtaining  the  midpoints  of  the  two  distorted  PN  sequences  with 
respect  to  the  master  word  group's  PN  sequences.  From  this 
computation,  the  slight  drift  between  the  samples  of  the  two 
tapes  can  be  calculated.  Using  this  drift  and  the  PN  sequence  » 
midpoints  of  the  master  tape,  an  estimate  of  the  midpoints  of  the 
PN  sequences  of  the  distorted  tape  can  be  made  taking  into 
account  the  shift  between  the  two  tapes.  The  true  midpoints  of 
the  PN  sequences  of  the  distorted  tape  with  respect  to  the  master 
tape  can  now  be  found  by  again  computing  the  cross-correlation 
between  each  pair  of  PN  sequences.  Using  the  distances  between 
the  midpoints  of  the  PN  sequences  and  the  words  following  them  of 
the  master  word  group,  and  the  drift  between  the  two  tapes,  the 
midpoints  of  the  words  of  the  distorted  word  group  with  respect 
to  the  master  can  be  obtained  using  the  midpoints  of  the  PN 
sequences  of  the  distorted  word  group.  This  procedure  is 
repeated  for  all  eight  word  groups  of  each  distorted  tape. 

The  alignment  of  two  words  from  two  different  tapes  to 
within  10  samples  was  the  goal  of  the  synchronization  procedure. 
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This  can  safely  be  assumed  using  the  above  alignment  scheme.  The 
PN  sequences  on  either  side  of  each  word  are  lined  up  to  within 
one  sample  in  all  cases.  The  drift  between  two  consecutive  PN 
sequences  was  never  more  than  15  samples,  usually  quite  a bit 
less.  Since  the  word  to  be  aligned  is  roughly  midway  between  the 
two  PN  sequences  and  the  drift  is  taken  into  account,  the  10 
sample  synchronization  specification  is  always  met,  generally  to 
within  a sample  or  two.  The  drift  between  two  tapes  was  verified 
to  be  linear,  with  only  small  (+1  sample)  fluctuations. 

Three  computer  programs  were  used  in  the  synchronization 
procedure  discussed  above.  "Words"  was  the  program  that  was  used 
to  find  the  locations  of  the  PN  sequences  and  words  through  the 
SD  energy  criterion.  FFTC0R4  computed  the  cross-correlation 
between  the  PN  sequences  of  the  two  tapes.  Finally,  WRDMIDP 
calculated  the  midpoints  of  the  words  of  the  distorted  tape  with 
respect  to  the  master  tape.  A listing  of  all  three  programs  can 
be  found  in  Appendix  A. 
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4.  DATA  PROCESSING  OF  THE  VOICE  TAPES 


Once  the  master  analog  tape  of  eight  50-word  groups  was  made, 
it  could  then  be  sent  over  various  voice  communication  channels 
to  obtain  distorted  tapes.  Copies  of  the  distorted  tapes  were 
sent  to  Fort  Huachuca  to  be  scored  for  intelligibility.  The 
intelligibility  score  for  a single  word  was  the  percentage  of  the 
listener  panel  that  correctly  identified  it,  and  the  intelligibil- 
ity score  for  the  entire  word  group  was  the  average  score  of  all 
50  words  in  the  word  group.  The  subjective  intelligibility 
scores  were  used  later  for  comparison  with  the  objective 
intelligibility  measure.  The  distorted  tapes  and  the  master  tape 
were  then  processed  to  obtain  the  LPC  information  necessary  to 
develop  the  objective  distance  measure.  No  filtering  was  used  on 
the  tapes  used  for  subjective  scoring. 
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A.  Digitization  of  the  Voice  Tapes 

Before  the  master  or  distorted  tapes  were  digitized,  they 
were  first  sent  through  a pre-emphasis  filter,  and  then  low- 
pass  filtered  to  3.2  kHz.  Pre-er.phasis  was  used  because  it 
enhances  the  high-frequency  formants  of  the  speech  which  is 
important  for  speech  comparison.  Pre-emphasis  also  limits  the 
effects  of  the  glottal  waveform  and  lip  radiation  and  therefore 
enhances  the  spectral  properties  of  the  speech  due  to  the  vocal 
tract. 

Based  upon  an  average  vocal  tract  length  of  17  cm,  the  first 
three  formant  frequences  will  be  found  in  the  frequency  range  of 
about  250  - 2800  Hz.  Shorter  vocal  tracts  will  shift  this  range 
up  slightly.  Low-pass  filtering  at  3.2  kHz  would  therefore  pass 
the  first  three  formant  frequencies.  This  (3.2  kHz)  is  generally 
also  the  high  frequency  cut-off  for  most  voice  communication 
channels.  Any  noise  above  3.2  kHz  picked  up  by  the  distorted 
tapes  will  also  be  filtered  out  which  will  help  the  accuracy  of 
the  objective  distance  measure. 

The  analog  tapes  were  sampled  at  10  kHz  and  then  quantized 
to  12  bits.  The  sampled  signal  was  then  stored  on  digital 
magnetic  tape  for  future  processing. 

B.  Analysis  Conditions 

Once  the  tapes  are  digitized  and  the  distorted  tapes  are 
aligned  to  the  master  tape,  LPC  processing  of  each  word  can  be 
done.  First,  however,  the  decision  must  be  made  regarding  which 
least  squares  method  to  be  used  to  obtain  the  predictor 
coefficients.  For  this  study,  the  autocorrelation  method  of 
linear  prediction  was  chosen  over  the  covariance  method  because 
it  requires  fewer  calculations,  it  is  assured  of  producing  a 
stable  filter,  (i.e.  all  the  poles  are  within  the  unit  circle), 
and  it  allows  a meaningful  spectral  matching  term  to  be  computed. 
Also,  as  mentioned  before,  the  autocorrelation  method  can  be 
considered  as  estimating  the  predictor  coefficients  by  approxi- 
mating the  spectrum  of  the  speech  waveform  by  an  all-pole  filter 


or  model.  A detailed  comparision  between  the  two  methods  is 
given  in  Makhoul  and  Wolf  (1972) . 

The  analysis  interval  in  which  the  speech  waveform's  spec- 
trum is  estimated  should  be  short  enough  so  that  vocal  tract 
movement  is  negligible,  but  long  enough  to  insure  stable  spectral 
estimates.  The  vocal  tract  can,  in  general,  be  assumed  to  be 
stationary  on  the  order  of  15  to  20  ms.  Since  in  the  auto- 
correlation method  of  linear  prediction  the  approximation  is  to 
model  a short-term  signal  spectrum,  it  is  necessary  to  window  in 
order  to  guarantee  spectrally  accurate  results.  By  using  a non- 
rectangular  window,  a larger  analysis  interval  can  be  used 
without  sacrificing  spectral  accuracy.  Therefore,  the  length  of 
the  analysis  frame  was  chosen  to  contain  256  samples,  which  means 
an  analysis  interval  of  25.6  ms  because  of  the  10  kHz  sampling 
rate.  The  non-rectangular  window  used  on  each  analysis  frcune  was 
a Hamming  window  of  the  form 

Wjj(129-n)=W^(128+n)=j-^jJ^.  (5.04+0.46  cos  ^)n=l,  2,  . . .,128 

(4.1) 

where  N=256,  the  analysis  frame  length.  The  Hamming  window 
was  chosen  because  of  its  desirable  spectral  properties  and  its 
widespread  use  in  linear  prediction  literature. 

Another  important  consideration  was  the  number  of  predictor 
coefficients  to  be  used.  As  a practical  matter,  it  is  best  to 
choose  the  niunber  of  coefficients  as  small  as  possible  because  it 
saves  computation  time  and  there  is  less  chance  of  filter 
instability  due  to  finite  arithmetic  effects.  It  has  been  shown 
that  to  represent  adequately  the  vocal  tract  under  ideal  circum- 
stances, the  memory  of  the  model  or  filter  must  be  equal  to' twice 
the  time  required  for  sound  waves  to  travel  from  the  glottis  to 
the  lips,  i.e.,  M = 2L/C,  where  L is  the  length  of  the  vocal 
tract  and  C is  the  speed  of  sound.  Using  the  average  vocal  tract 
length  of  L =»  17  cm  and  the  speed  of  sound,  C = 34  cm/ms,  the 
memory  required  is  1 ms  with  a 10  kHz  sampling  rate,  the  number 
of  predictor  coefficients  needed  is  equal  to  the  sampling  rate 
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times  the  memory  required,  or  10  coefficients.  In  order  to  take 
into  account  the  influences  of  the  glottal  waveform  flow  and  lip 
radiation  characteristics  an  additional  two  coefficients  are 
necessary.  Hence,  twelve  predictor  coefficients  were  computed 
for  each  analysis  frame  thereby  producing  a twelve  pole  filter 
that  accurately  models  the  spectral  properties  of  each  speech 
interval. 

The  information  necessary  to  choose  the  various  parauneters 
and  methods  for  LPC  processing  were  obtained  from  experimental 
results  and  several  references,  which  include  Markel  and  Gray, 
(1973)  and  Makhoul  (1973),  and  Boll  (1973). 

C.  LPC  Data 

With  all  the  necessary  parauneters  and  methods  chosen  to 
analyze  the  speech,  LPC  processing  of  each  word  could  then  be 
undertaken.  Computer  program  LPC  was  used  to  perform  the  speech 
analysis.  A listing  of  the  program  can  be  found  in  Appendix  A. 
The  256  point  window  was  moved  along  each  word  at  256  sample 
shifts.  Originally,  the  window  was  shifted  by  128  points  to 
create  overlapping  analysis  fraunes.  However,  due  to  the 
averaging,  significant  differences  in  the  distance  measures,  to 
be  described  in  Chapter  5 .were  not  found  when  the  two  methods 
were  compared.  Consequently,  the  256  point  shift  was  adopted  to 
save  processing  time.  The  distorted  word  lengths  were  naturally 
defined  to  be  the  same  as  the  corresponding  master  word  lengths. 

Processing  of  the  windowed  signal  to  obtain  the  LPC  para- 
meters was  done  using  the  Levinson  Algorithm.  A flow  chart  of 
the  algorithm  is  shown  in  Figure  4-1.  In  order  to  start  the 
algorithm,  the  normalized  autocorrelation  terms  of  the  windowed 
signal  had  to  be  obtained.  This  was  accomplished  using  the 
direct  method  of  equation  (2.7)  as  opposed  to  using  the  Fast 
Fourier  Transform  (FFT)  because  of  the  small  number  of  lag  terms 
needed  (thirteen  terms) . The  Levinson  Algorithm  was  used  to 
obtain  the  twelve  predictor  coefficients,  (a^^)  , and  the  minimum 
total  squared  error,  Each  reflection  coefficient  (k^)  was 
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Figure  4-1.  Flow  diagram  for  LPC  processing 


checked  to  insure  that  it  was  less  than  0.99/  thereby  guarantee- 
ing a stable  filter. 

All  of  the  LPC  information  was  then  stored  on  magnetic  tape 
to  be  used  later  in  comparing  the  speech.  The  LPC  data  stored 
consisted  of  R^,  the  autocorrelation  normalization  factor  and 
energy  term  for  the  analysis  frame,  the  normalized  autocorrelation 
terms,  (r^^)  , the  predictor  coefficients,  (a^^)  , the  reflection 
coefficients,  (k^)  , the  total  minimum  squared  error,  and 

the  autocorrelation  terms  of  the  predictor  coefficients,  (g^^)  . 

The  (g^)  terms  simplify  the  calculation  of  the  distance  measures 
to  be  discussed  in  the  next  chapter. 

5.  DISTANCE  MEASURES 


The  distance  or  distortion  measure  to  be  developed  is  inten- 
ded to  predict  accurately  differences  between  two  intervals  of 
speech.  One  of  the  speech  intervals  will  always  be  from  the 
undistorted  speech,  while  the  other  will  be  a distorted  version 
of  the  original  speech.  The  distance  measure  obtained  must 
therefore  be  interrupted  relative  to  distance  measures  of  other 
distorted  speech  segments  when  compared  to  the  undistorted 
version.  The  distance  measure  is  then  mapped  into  an  objective 
determination  of  the  intelligibility  of  the  speech.  The  distance 
measures  discussed  in  this  chapter  are  all  developed  in  the 
framework  of  LPC. 


A.  Linear  Prediction  Residual  Distance  Measures 

To  digress  a moment,  a sample  is  estimated  by  the  linear 

combination  of  the  preceeding  "P"  samples  in  LPC.  The  error  or 

linear  prediction  residual  can  then  be  written  as 

P 


e 


n 


S 


n 


E 


a . 
1. 


S 


n-1  ' 


(5.1) 


where  (a^^)  are  the  predictor  coefficients.  These  coefficients 
are  obtained  by  choosing  them  so  as  to  minimize  the  total  squared 
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error.  The  total  squared  error  or  linear  prediction  residual 
energy  can  be  considered  to  be  the  output  of  an  inverse  filter 
H(z)^,  where 

i=l 

[H(z)]^  is  the  filter  that  minimizes  the  residual  energy  and  H(z) 
corresponds  to  a smoothed  spectral  estimate  of  the  data  sequence 
(Sy^)  up  to  a scale  factor  representing  the  gain. 

If  (S^)  is  passed  through  a different  inverse  filter, 
iH'(z)]  of  the  form 


[h'(z)]"^  = 1 + 


E alz- 

i=l 


(5.3) 


which  minimizes  the  residual  energy  for  some  other  data  sequence 

I 

(S^  ) , the  residual  energy  D,  must  be  greater  than  or  equal  to 
the  minimum  residual  energy  E,  i.e.,  D^E,  with  the  equality  hold- 
ing  if  and  only  if  H(z)  = H (z) . Assuming  the  data  sequence  (S^) 
is  obtained  from  an  analysis  frame  of  speech  from  the  undistorted 

t 

tape  and  (S^  ) is  the  corresponding  analysis  frame  from  a dis- 
torted tape  the  difference  between  D and  E is  a measure  of  the 
distance  between  the  two  speech  segments.  (Unless  otherwise 
identified,  unprimed  variables  represent  data  from  the  master  or 
undistorted  tape.) 

The  dual  of  the  above  situation  is  also  true.  If  (S  ) is 

1 I I “ 

sent  through  tH(z)]  , the  output  will  be  D , while  E is  the 

• 1 • • • 

output  of  [H  (z)]  when  (S_  ) is  sent  through  it.  Again,  D >E 

n I 

with  equality  if  and  only  if  H (z)  = H(z).  As  before,  the 

I I 

difference  between  D and  E can  be  considered  to  be  a distance 
measure  between  the  two  speech  segments. 

I t 

E,  E , D,  and  D can  all  be  written  as  a combination  of  the 
autocorrelation  terms  of  (S^)  and  (S„ ' ) and  the  corresponding 
linear  prediction  coefficients  (a^)  and  (a^) . Let 


h.  ~ ~®i'  ~^2'  * ' 


(5.4) 


be  the  transpose  of  the  linear  prediction  coefficient  vector  A, 
and 


R = n . 


i-k 


i/k  — 0^  1/  • • m f P 


(5.5) 


the  normalized  autocorrelation  matrix.  The  four  error  terms  can 
then  be  written  as 


E = a"^  R A 


(5.6) 


■ Im  I I 

E = A R A 


I IP  I 

D = A R A 


T • 

' = A R A 


(5.7) 

(5.8) 

(5.9) 


where  the  primes  signify  variables  from  a distorted  tape.  The 
derivation  of  this  can  be  found  in  Market  and  Gray  (1973)  and 
Boll  (1974). 

f 

E and  E are  calculated  for  each  analysis  frame  through  the 

I 

Levinson  Algorithm,  but  D and  D are  calculated  when  a distorted 
tape  is  compared  to  the  master  tape.  These  calculations  can  be 

I 

simplified  because  of  the  structure  of  R and  R_,  the  auto- 
correlation matrices,  by  calculating  the  autocorrelation  terms  of 

f 

A and  A , the  linear  prediction  coefficient  vectors.  Using  the 

f I 

symmetry  of  the  autocorrelation  terms  of  A,  and  D and  D , can 
be  written  as 

- t 


(5.10) 


i=0 
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(5.11) 


where 


i=0 

• P-lil 

= 2 . 2 ajs+i 

k=0 


(i  = 1,  2,  . . . , P 


(5.12) 
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k=0 


(5.13) 


(5.14) 


(5.15) 


(a.)  and  (a^)  are  the  P+1  terms  in  the  vectors  A and  ^ 
respectively.  (9j_)  and  (gj_)  can  be  calculated  for  each  analysis 
frame  right  after  the  predictor  coefficients  are  obtained  and 
stored  on  the  LPC  data  tape. 

Each  of  the  four  error  terms  can  be  interpreted  in  the 
frequency  domain.  Using  Parseval's  Theorem,  the  total  squared 
error,  E,  can  be  written  as 

.tt/T 


= S en2  = ^ j 


|E(a))  I do)  , 


ib)t 


(5.16) 


where  E(u))  is  obtained  by  substituting  z * e into  E(z).  From 
Chapter  2,  the  minimum  linear  prediction  error  was  found  to  be 


E(z).  = S(z)  lH(z)  ]“^ 

while  the  least  squares  estimate  can  be  written  as 


(5.17) 


E(z)  = S(z)  [H(z)  ] 


-1 


(5.18) 


substituting  z = into  (5.17)  and  (5.18),  one  obtains 

E(u))  « 8(0))  [H(o))]"^  (5.19) 
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(5.20) 


^ /V 

E(a))  = S(u)  [H(a))  . 


Rearranging  (5.20)  and  substituting  in  E(z)  = A/ 

[H(a))]-^  = . 

S(w) 

Inserting  (5.21)  into  (5.19), 


(5.21) 


E(cl))  = A 

3(0)) 

substituting  (5.22)  into  (5.16) 


(5.22) 


|S(a)) 


(5.23) 


2 ^2 

But  |S(o))l  and  (S(o))|  are  just  the  corresponding  power  spectra, 
P(o!)  and  P(o)),  of  the  speech  signal  and  its  least  squares  linear 
prediction  estimate.  Therefore, 


do) 


(5.24) 


P(o)) 


Similarly  E can  be  shown  to  be 


' 2 

p’  ta 

-tt/T  , 

1 P (O)) 

^ ■ -TT 

can  also  be 

obtained  by 
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D - 

-tt/t 

1 P (O)) 

° 27r 

L/t 

' TA^ 

° = IT 

-tt/T  , 

j p (O)) 

J-»/T 

(5.25) 


(5.26) 


(5.27) 


The  distance  measures  D and  D are  not  all  that  pleasing 
when  defined  in  the  frequency  domain.  They  compare  the  ratio 
differences  between  a true  speech  power  spectrum  and  estimated 
power  spectrum.  A much  more  desirable  measure  would  compare  the 
ratio  differences  between  the  estimated  power  spectra  of  the 
undistorted  speech  and  the  distorted  version  of  it.  This  can  be 

t I 

done  by  taking  the  ratio  of  D to  E and  D to  E . The  ratio  of 

f I 

each  of  these  pairs  of  residual  errors,  D/E  and  D /E  , then 
defines  two  new  distance  measures  which  are  much  more  appropriate. 
In  both  cases,  the  ratios  are  greater  than  or  equal  to  one,  with 

I 

equality  if  and  only  if  H(z)  = H (z) . 

f I 

The  ratios  D/E  and  D /E  are  sometimes  called  likelihood 
ratios  because  under  certain  circumstances,  they  have  been  shown 
to  be  true  likelihood  ratios  by  Itakura  (1975) . As  mentioned 
before,  the  frequency  domain  interpretation  of  the  likelihood 
ratios  gives  a good  justification  for  using  them  as  distance 
measures.  In  the  time  domain 


(5.28) 


where 


E 


(5.29) 

(5.30) 
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Gray  and  Markel 

(1976)  have  shown 

that  D/E 

frequency  domain 
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where  the  substitution  z = is  made  in  the  filters  H{z)  and 

I 

H (z) . Inverting  (5.21)  and  its  dual,  one  obtains 


H(a))  = 


- 5(0)) 


h'  (o))  = g- 


(5.32) 


(5.33) 


Substituting  (5.32)  and  (5.33)  into  (5.31) 

.A 


D _ T 
E JtT 


S(o)) 


-n/T 


Is'  (0))  1^ 


(5.34) 


f 

l./T  Is  ««) 


(5.35) 


Once  again,  the  magnitude  squared  of  the  signal's  spectrum  is 
just  its  power  spectrum,  therefore 


•'-i 


Similarly 


P ((ij) 


(5.36) 


P (o)) 


(5.37) 


-n/T 


As  can  be  seen  from  (5.36)  and  (5.37),  D/E  and  D /E  compute  the 
differences  between  the  estimate  power  spectra  of  the  undistorted 
and  distorted  speech,  while  D and  D compared  the  true  power 
spectra  of  the  speech  to  estimates  of  the  power  spectra. 

I I 

D/E  and  D /E  will  be  used  as  the  basic  components  in  the 

distance  measures  discussed  in  B,  below.  Several  methods  are  used 

« 1 

to  normalize  D/E  and  D /E  so  that  two  distorted  tapes  can  be 
compared  relative  to  each  other  by  comparing  them  only,  to  the 
master  or  undistorted  tape.  For  more  information  concerning  the 
lilcelihood  ratios  above,  see  Gray  and  MarJcel  (1976)  . 


One  additional  measure  that  is  considered  here  is  derived 

Then, 


as  follows.  Let  S + N = S„ 

n n N 


5 -Ln^n-i 


(5.38) 


This  can  be  rewritten  as 

I I 

E + + [cross  products]  = D , 

I ip 

where  is  given  by  A with  the  ^ the  autocorrelation  matrix 
of  the  noise.  Assuming  the  signal  and  the  noise  are 
uncorrelated,  this  simplifies  to 


E+D„'.d’.  ■ (5.39) 

Assuming  the  noise  for  any  frame  (N^)  is  the  same  (statistically) 
except  for  a constant  factor  as  the  noise  for  a period  during 
which  voice  is  not  present  (N^) , we  have 


I 


D 


N 


(5.40) 


where  k is  a different  constant  for  each  frame.  Since  k must  be 
positive,  we  have 

^ ^ Id'  - El  (5.41) 


A large  number  of  calculations  showed  that  E/D  „ was  very  small 

• I ” 

compared  to  D /D  jj,  and  consequently,  k was  calculated  using 

k = d'/d'~  (5.42) 

k was  then  averaged  over  all  freimes  of  all  words  to  obtain  Ic. 
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A signal  to  noise  ratio  was  determined  using  the  peak  signals 
during  a PN  sequence  and  the  noise  N from  a quiet  period  between 
PN  sequences,  for  each  word  group.  The  quantity 


SNR  = 10  logj^Q  S/N  - 10  logj^Q  It 


(5.43) 


was  then  used  to  calculate  an  AI  score, 
0 SNR< 0 

AI  = J if  0<SNR<30 

^1  SNR>30. 


(5.44) 


B.  An  Objective  Intelligibility  Measure 

I I 

The  quantities  D /E  and  D/E  were  computed  on  a frame-by- 
frame  basis  using  the  computer  program  DISTMEA  and  stored  along 
with  other  LPC  data  for  use  in  developing  an  objective  intelligi- 
bility measure.  (See  Appendix  A.)  The  natural  logarithms  of 

• I 

D /E  and  D/E  are  respectively  labeled  El  and  E2  and  relate 
directly  to  a decibel  (dB)  scale. 

I 

Under  the  assumption  that  the  errors  e ^ are  independent 

Gaussian  variables,  Itakura  derives  the  result  that  El  is  a 

chi-squared  variable  with  P(=12)  degrees  of  freedom.  Here, 

because  of  the  windowing,  101.  However,  since  in  general 

the  e are  correlated,  it  is  assumed  that  the  actual  N is 
n eft 

smaller  than  this  (101) . 

Instead  of  modifying  the  procedure  given  in  the  next 

paragraphs  was  used  to  modify  two  thresholds.  First,  a lower 
threshold  for  El  was  taken  as  0.82,  based  on  an  average  of 
approximately  three  times  the  "barely  perceptible"  difference  of 
Flanagan  (1972 ) and  three  times  the  "barely  perceptible"  thres- 
hold used  by  Scimbur  and  Jayant  (1976)  . Values  of  El  below 
0.82  mean  the  frame  is  understood.  An  upper  threshold  of  2.46 
(3x0.82)  was  used  to  decide  that  the  frame  was  completely  mis- 
understood. A linear  relationship  was  used  between  0.82  and 
2.46. 


32 


i 


A scunple  of  noise  was  taken  from  the  distorted  tape  being 
analyzed  (actually  the  same  samples  used  in  the  previous  section 
to  derive  k) . From  this  sample  several  frames  were  analyzed  from 
which  two  frames  were  selected  using  the  criteria  of  the  largest 
and  smallest  values  for  the  sum  of  the  squares  of  the  predictor 
coefficients,  these  two  frames  were  used  with  the  master  tape  to 
calculate  (frame  by  frame)  values  of  EIN^  and  E1N2»  where  N 
signifies  noise.  These  two  values  (for  each  frame)  were  then 
averaged  to  obtain  EIN.  If  EIN  <2.46  the  thresholds  were  not 
changed.  If  EIN  >2.46  the  thresholds  were  changed  to 

0.82  + 0.82  (EIN  - 2.46)  and 
2.46  + 0.82  (EIN  - 2.46)  . 

To  summarize,  two  thresholds  T1  and  T2  are  defined.  Using  these, 
a linear  measure  is  defined  for  each  frame  as 


LMl  = 1 if  El  < T1 

= 0 if  El  > T2 


(5.45) 


T2-E1 

T2-T1 


otherwise . 


The  linear  measure  LMl  of  the  above  method  was  calculated 

for  all  frames  of  each  word.  Further,  for  each  word  an  average 

of  this  measure  was  calculated  for  those  frames  for  which 

Rq^Rq/2,  where  R^^  was  the  average  of  R^  for  the  word.  This  was 

designated  LMIH.  Similarly,  an  average  for  frames  for  which 

R <R  /2  was  calculated  and  called  LMIL.  This  divides  the  measure 
o o 

into  two  groups,  one  for  frames  with  higher  power  and  one  for 
frames  with  lower  power  compared  to  the  average  power  in  the 
word.  Frequently,  although  not  always,  the  low  power  frames 
correspond  to  the  unvoiced  speech  and  the  high  power  to  the 
voiced  speech.  Finally,  LMIH  and  LMIL  were  averaged  over  fifty 
words  to  form  LMIH  and  EHTE. 

In  a similar  fashion,  high  and  low  values  EIH  and  ETX, 
and  eInh  and  EINL  were  calculated  in  order  to  modify  the 


M 


Va?**  ? 
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average  of  the  linear  measures  in  the  following  way.  If  EIH 
(resp  eIl)  is  less  than  .82  no  modification  is  made.  If  EINH 
(resp  EINL)  is  greater  than  2.46  no  modification  is  made. 
Otherwise  LMIH  (resp  lMIL)  is  multiplied  by 


EINH  - 0.82 

2.46  - 0.82 


(resp 


EINL  - 0.82 
2.46  - 0.82 


) . 


The  modified  high  and  low  measures  were  then  averaged  to 
form  LMl.  This  has  the  effect  of  weighting  the  low  values  more 
than  the  high  values  since  only  about  ]/3  of  the  frames  are  low. 

I 

The  correlation  Cj^  , of  and  was  then  calculated, 
and  multiplied  by  LMl  t8  ?orm  CLMl.  This  value  was  averaged  with 
the  AI  measure  of  the  previous  paragraphs  to  form  ASQ,  the 
objective  articulation  score. 

The  correlation  correction  was  applied  to  account  for  fading 
signals.  Several  other  methods  were  used  which  compared  the 
signal  levels  on  a word-by-word  and  frame-by-frame  basis.  These 
required  considerably  more  computing  time  and  gave  essentially 
the  same  results  as  using  the  correlation  factor. 

The  quantity  AI  is  shown  in  Figure  5-1  plotted  vs  the 
subjective  articulation  scores  (AS) , the  quantity  CLMl  is  shown 
in  Figure  5-2  plotted  vs  AS,  and  ASO,  the  average  is  shown  in 
Figure  5-3  plotted  vs  AS.  In  Figure  5-3,  the  bars  indicate  the 
confidence  limits  about  the  subjective  score. 

While  subjective  intelligibility  scores  for  isolated  words 
are  repeatable,  there  are  some  fluctuations  in  the  scores.  A 
listener  panel  is  trained  on  a set  of  word  groups  with  well 
established  intelligibility  scores  and  standard  deviations.  The 
average  intelligibility  of  each  training  word  group  scored  by 
individuals  that  make  up  a listener  panel  is  always  plus  or  minus 
one  standard  deviation  of  its  actual  intelligibility  score. 

Also,  the  standard  deviation  of  the  listener  panel  for  each 
training  word  group  is  approximately  the  seune  as  its  actual 
standard  deviation.  Using  the  training  procedure,  a listener 
panel  will  produce  an  intelligibility  score  within  plus  or  minus 
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two  standard  deviations  of  the  actual  intelligibility  score  95% 
of  the  time. 

CLMl'  was  calculated  for  7 additional  tapes  (each  containing 
eight  50-word  groups) . For  these  calculations,  a different  noise 
frame  was  used  (see  Section  5 A ) for  each  word.  A noise  frame 
was  chosen  immediately  following  the  PN  sequence  before  the  word 
being  evaluated.  These  frames  were  used  to  calculate  EIN  for  the 
words.  These  data  are  shown  in  Figure  5-4.  Two  points, 
represented  by  circles  are  suspect.  The  analog  tapes  used  for 
the  data  analyses  were  transcribed  onto  1/4  inch  tapes  for 
subjective  scoring.  Two  word  groups  were  somehow  omitted  from 
the  1/4  inch  tape  for  this  system  and  two  other  word  groups 
(represented  by  the  circles)  gave  the  anomalous  scores.  Time  did 
not  permit  making  a new  tape  and  subjectively  scoring  it  to 
determine  the  validity  of  these  scores.  The  subjective  scores 
for  the  entire  tape  are  shown  in  Table  5-1. 

The  objective  measures  CLMI  shown  in  Figures  5-2  and  5-4  are 
essentially  the  same  calculation,  varying  only  in  the  detail  of 
whether  one  or  many  noise  frames  were  used  i i the  calculations. 
The  combination  of  all  the  points  (90)  in  th 2se  two  figures, 
excluding  the  two  circled  points  yields  a cc rrelation  of  .982 
between  CLMI  and  AS.  The  RMS  error  between  AS  and  CLMI  is  5.5(%). 


Table  5-1 
Tape  T3 

Word  Group  AS  % 


361 

60.70 

312 

68.00 

291 

70.50 

265 

72.70 

275 

No  score 

305 

• 

No  score 

214 

45.00 

283 

17.20 

Figure  5-4.  A comparison  of  the  objective  measure  CLMl, 
calculated  using  a noise  term  for  each  word 
with  the  subjective  measure  AS. 


6.  DISCUSSION  AND  CONCLUSIONS 


i 

I 


i 
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Based  on  the  results  of  the  feasibility  study  (Hartman  and 
Boll,  1976)  several  refinements  and  normalizations  of  the  distance 
measures  have  been  ex2unined.  Throughout  this  study,  the  under- 
lying philosophy  has  been  to  choose  measures  requiring  as  few 
calculations  as  possible  while  still  producing  a reasonable 
correlation  with  the  subjective  scores.  Several  of  these  measures 
are  displayed  here,  namely  CLMl  using  the  same  noise  term  for  all 
words  on  a tape,  CLMl  using  a different  noise  term  for  each  word, 
and  AI.  The  best  correlation  is  obtained  using  CLMl  calculated 
using  the  noise  term  close  to  the  word  as  might  be  expected, 
particularly  if  the  noise  is  changing  in  character.  Since  using 
the  different  noise  terms  does  not  significantly  increase  the 
number  of  calculations  required,  this  CLMl  measure  appears  to  be 
the  prime  candidate  for  hardware  implementation. 

In  the  software  implementation,  about  70%  of  the  required 
computer  time  is  used  in  the  alignment  procedures,  which  also 
accounts  for  about  85%  of  the  manpower  requirements.  Of  the 
remaining  30%  of  computer  time,  about  70%  (21%  of  the  total)  is 
required  for  obtaining  the  LPC  parameters  and  the  rest  (9%  of 
the  total)  is  required  for  computing  the  distance  measures  which 
are  displayed  in  this  report. 

The  synchronization  developed  here  is  adequate  for  word 
alignment  even  when  severe  distortion  is  introduced  by  the  system 
being  tested.  With  a view  to  hardware  simplification,  the  PN 
cross  correlations  were  computed  using  one-bit  quantizing  and  the 
alignment  was  within  +1  Scunple  of  the  result  obtained  using  the 
12  bit  quantizing.  This  means  that  binary  devices  can  be  used  in 
a hardware  synchronization  circuit.  Consequently,  since  real 
time  LPC  devices  have  been  built,  the  feasibility  of  a hardware 
implementation  is  established. 

Although  not  all  types  of  voice  communications  systems  were 
studied  in  this  effort,  representative  voice  systems  were  used. 

The  results  of  the  study  indicate  that  the  objective  measure (s) 
developed  are  good  predictors  of  the  subjective  scores. 
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In  order  to  enlarge  the  data  base  further,  or  to  investigate 
other  uses  and  modifications  of  these  methods  it  appears  that  the 
most  economical  procedure  is  to  develop  a flexible  hardware  system. 
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APPENDIX 


This  appendix  contains  the  program  lists  for  processing  the 
voice  data.  Table  A-1  shows  the  flow  diagram  for  the  use  of  the 
programs. 

The  lists  are  arranged  with  all  of  the  progreuns  first, 
followed  by  the  subroutines  used  in  the  programs.  Each  program  or 
subroutine  is  identified  by  name  and  purpose  with  ccxnment  state- 
ments at  the  beginning. 

Commonly  used  library  subroutines  are  not  included  in  the 
subroutine  listings,  and  some  of  these  may  need  slight  revision 
when  used  with  a different  computer  system.  Particular  attention 
should  be  given  the  input-output  statements  and  instructions 
dealing  with  packing  data. 
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BRO&RAM  ilQRDS(IMPUT<OUTPUT«T«Pei> 

PROGRRN  WORDS  PINOS  THE  LOCATION  OP  THE  PN  SEIUENCES  AND  WOTOS  OP 
A BARTICULAR  HOID  GROJP  PY  FINDING  THE  HEAN  ANO  STANDARD  OEVIATION 
OF  CONSECUTIVE  125  SAHPLL  RLOCRS.  THE  STANOARO  DEVIATION  IS  USED 
AS  AN  ENERGY  CRITERION  IN  ORDER  TO  LOCATE  THE  PN  SE3UENCES  AND  WORDS. 
OINENSION  INISOl  tlRECOROIAOO>.IO(lP5),ITEH'»<2J#0) 

initialize  variables.  •THOLOt  IS  A THRCSHOLO  USED  TO  OETERNINE  IP 
BART  OF  A PN  SE2UENCE  OR  WORD  IS  PRESENT  IN  THE  DATA.  .ISKIP*  AND 
•lENO*  OETERHINE  THE  STARTING  AND  ENDING  »OINTS  0«  THE  RECORDS  TO 
BE  “ROCESSEO.  •ROUNT.,  aICYCLE*.  AN"  .IPARTS*  ARE  COUNTING  VARIABLES. 
»ILENGTH»  is  the  NUHBER  op  SANBLES  TO  BE  “RDCiSSED  AT  A Tiwr  twD 
•IWOROS*  IS  THE  NUHBER  OF  SAHPLES  IN  A TA»E  RECORD. 

TH0LD=32.0 

IS<IP=r 

IEN0=3R5' 

KOUNT=l 

ICYCLEsU 

I»ARTSsl6 

ILENGTH=125 

IWOROS=A  .0 

IREC^ISHIP 

the  tape  is  positioned  TO  THE  STARTING  POINT  FOR  PROCESSING. 

IF  (ISKIP  .CD.  0)  60  TO  E 

DO  1 RIOH.ISKIP 

»JFFER  IN  <1»1)  IINID.INCB^)) 

IF  lUNITIDt  lOtCOtlR 
1C  CONTINUE 
5 BRINT  55»ISK1» 

55  bqrhaT  IlHl.IS.*  RECORDS  SKIP»ED  •) 

I'  ?REC=IREC*3 

FIVE  CONSECUTIVE  Sr  WORD  RECORDS  ART  BUFFERED  IN  ^OR  PROCESSING. 

EACH  6C  BIT  WORD  IS  UNBACKED  INTO  =IVE  12  BIT  WORDS  USING  SJBRnuIINE 
JNBACK.  THIS  YlILDS  A TOTAL  OP  2000  SAHPUES  TO  HE  PROCESSED  AT  A TIME. 
BJPPER  IN  <1,1»  «INIl>tIN<B  )) 

IP  lUNITID)  I5,20»20 
25  brinT  12C.IREC 

J2D  PORNAT  I.  RECORD  ..I5.*  HAS  A “ARTTY  ERR'R  •) 

GO  TO  3 

15  LLiLENGWd) 

CALL  UNPACK  CIRECORD.IN.LLl 
KINDEX=ICYrLE*IWDHDS 
DO  1'  K12il,lW0RDS 
12  ITE'”<K12«KINDE«)sIREC0RD»k;2> 

ICYCL£=ICYCLE»1 

I'  IICYCLE  .LT.  G)  SO  TO  .'I 

THE  2000  SAHolES  ARE  3R1wrN  UP  INTO  SIYTeEN  I’S  SAMPLE  ^LDCKS.  POo 
EACH  3LDCK.  THE  “EAN  AND  STANDARD  DEVTATIDN  /RE  CDHpUTED.  IP  THE 
STANCARD  DEVIATION  IS  GREATER  THAN  OB  EDUIL  ID  THE  PRESET  THRESHOLD. 
THOLD.  THE  FIRST  TEN  SAHPLES.  Tnr  <>.;an.  THE  STANDARD  DEVIATION.  AND 
The  “LOCK  NUH3ER  pOS  TH"  SAMPLE  “LOCK  ARF  PRINTED.  OTHERWISE  THE 
PR35RAH  PROCEEDS  TO  THE  NEXT  SAHpuE  BLOCK. 

DO  35  K55  = l. IPARTS 
NiNj/AUKSs-:!.  I Length 
DO  *f  KAOsl.ILENGTH 
AO  lD<KA,)sITENPHKAC.NINOEX) 

SUMlsSUN2sO.? 

00  A5  KASsl.ILENGTH 
SUBlsSUHl*ID(RA5> 

AS  SUN2:SUN2*101RA5)*I0(KAS> 

AbEANbSUNIZILENGTH 

STDOEVeSORTIABSI «1LENGTH.SUH2-SUH1.SUN1)/ILENSTH..2> ) 

IF  ISTDOEV  .LT.  THOLOI  GO  TO  50 

PRINT  ITO.KOUNT.  (I0(  I >.I:1.1(iI.AHEAN.STD0EV.K0UNT 
130  PORNAT  (•  ••IG,ini9.2F12.A«lB> 

50  ROUNTaROUNT*! 

35  CONTINUE 
ICYCLEal 

The  NUHBER  OF  THE  NEXT  RECORD  TO  BE  PROCESSED  IS  COMPARED  TO  THE 
NUMBER  OF  THE  LAST  RECORD  TO  BE  PROCESSED.  lEVD.  IF  IT  IS  LESS  THAN 
lEVO.  THE  PROGRAH  CONTINUES.  OTHERWISE  THE  PROGRAM  TERMINATES. 

IF  IIREC  .LT.  TEND)  SO  To  3J 
20  PRINT  IGO 
100  format  I1H1> 

PRINT  115.IREC 


115  FORMATIISt 
BRINT  110 

110  fdbnaT  <•  NORMAL  termination  •> 
CNO 
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<>«06RAN  FFTC0R4  ( INPUTtOUTPUTtPUNCHt TA»C  1 • Tt PiZ  > 

C PA06RAN  FFTC04A  COMMUTES  THE  CR 3SS-C0RRELA TI OM  BETWEEN  TWO  REOUCNCCS 
C OF  OATA  samples*  which  in  this  case  are  the  >N  SraUENCES. 

COHHON  /FFTl  /A  lAZROI*  B IA200)«  C (A200>«  D (A?00>t  01  (5900) 

5 COHHON  /FFT2  /ICOR  (21001 

DIMENSION  10(503  I •ITEMP(lj0>«IN(100> 

C INITIALIZE  variables.  •MPRINT*  IS  A PRINTING  CONTROL  VARIABLE  ANO 

C >NUHREC*  IS  THE  NUMBER  OF  SAMPLES  IN  A TAP-  RECORO.  *1ASIZE1>  AND 

C •IASIZE2*  ARE  COUNTING  VARIABLES.  .NRECREO*  ANO  (NRECIAA.  NEEP 

10  C TRACK  OF  THE  POSITION  OF  THE  UNOISTOPTEO  ANO  JISTORTED  TAPES. 

MPRINTsO 
NUMRECsAOO 
lASIZEI  = A2C0 
IASIZE2  = 5CC0 

15  PRINT  150R 

150A  FORMAT  (IHll 
ICO  NRECREO  = NRECIAA  s 0 
135  DO  110  KIO  = 1.  lASIZEl 

A (KlOl  = B (KlOl  = C (Kini  = 0 (KlOl  = C. 

20  110  CONTINUE 

00  115  K20  = 1.  1ASIZE2 

01  (K20>  = 0 

115  CONTINUE 

C THE  CONTROL  CARO  IS  READ  THAT  GIVES  THE  H1D»0INTS  OF  THE  TWO  SEOUENCES 

25  C TO  BE  PROCESSED  ANO  THEIR  LENGTH.  IF  VARIABLE  NIDPNTI  IS  EOUAL  TO 

C ZERO.  THE  program  TERMINATES.  INDEACS  ARE  THEN  COMPUTED  THAT 

C POSITION  THE  TWO  TAPES  IN  ORDER  TO  GET  THE  OESIREO  DATA. 

RE  AO  1532.MIDPNT1.MIDPNT2.NUMEAO 
150?  F0»HAT(31M> 

SO  IF  (M?D»NT1  .EO.  0)  GD  TO  AV5 

NUMBER  = 2 . » NUMEA= 

INOEXl  = HIOPNTl  - NUMBER  / 2 

INDEX’  = MIDPNT’  - NUMBER  / 2 

ISTRTRl  = INDEXl  / NUMREC  - NRECREO 

35  ISTRTR?  = INDEX?  / NUMREC  - NRECIAA 

lORIGIN  = INDEXl  - NUMREC  • (INOEXl  / NUMREC) 

IDRIIAA  r INDEX?  - NUMREC  • (INDEX2,  ! NUMREC) 

C THE  TWO  TAPES  ARr.  OOSITICINEO  TO  THE  START  df  THE  SEQUENCES  TO  BE 

C PROCESSED. 

♦ 0 IFdSTRTRl  .LT.  I)  SO  TO  5^0 

DO  lOG  HAD  = It  ISTRTR: 
buffer  in  fl,l)  (ITEMP(t).ITFMP(RD)) 

IF  (JNIT(  1))1A),  AB5.  IBS 
133  PRINT  13)8.  KAO 

A3  1308  formst  I.RFAD  ERROR  TAPF  1 nC) 

lAD  CONTINUE 
ROD  CONTINUE 

NRECREO  F NRECRiD  ♦ ISTRTRl 
IFilSTRTRZ  .LT.  1)  SO  TD  305 
SC  DO  133  HG3  r I,  TST0TR2 

HUFF1.R  I)  (’.I)  (ITEMPdl.lTEMPIOi)) 

IF  (UNIT!  D))133.  A.A,  IE.’ 

IbC  PRINT  1310.  KGD 

1510  FORMAT  (.REAO  ERROR  TAPE  J *,  HD 
55  155  CONTINUE 

595  CONTINUE 

NRECIAA  = NRECIAA  ♦ ISTRTR2 

C THE  TWO  OATA  SEJUENCES  ARE  BUFFERED  INTO  THE  PROGRAM.  UNPACKED  BY 

C SUBROUTINE  UNPACK.  AND  STORED  IN  INDIVIDUAL  ARRAYS. 

GO  NRECIN  : 2 * NUMBER  ! NUMREC 

IOFFSET  f 0 

00  1R9  KlOO  F 1.  NRECIN 
buffer  in  (1.1)  (IN(lltIN(RC>> 

IF  (UNIT(  imro.  A85*  165 
65  165  KPARITT  = KAO  * KlOO  - 1 

PRINT  1516. KPARITT 

1516  FORMAT  (.PARITT  ERROR  IN  OATA  RECORD.  RUN  IBORTEO  *.  Ilfj) 

50  TO  ATS 
170  CONTINUE 

TO  LLfLENGTHH) 


CALL  UNPACKdO.lN.NUMREC.LL) 
00  175  KBS  = 1.  NUMREC 
D1  (K95  ♦ IOFFSET)  f IO  (KB3) 


B (<152) 
« (K1S2> 
B IK152> 


A (K152> 
B fK152l 


CONTINUE 

NAECAtO  : NRECSEO  ♦ NRECIN 

00  1«5  K102  = It  NUN“EA 

A *<in2»  = 01  (lOAIGIN  ♦ <n2l 
CONTINUE 

DO  no  AIDS  = It  lAOlZE.'^ 

01  <A1P!1)  = ■’ 

continue 

IO'ESET  = 

DO  :35  <1JC  : It  NRECIN 

IN  (2tl)  (INI  l)tlN(“ '.)) 

IT  (JNIT<  2»)2!5t  ABBt  205 
AOAAITr  T 560  ♦ KlOO  - 1 
PAINT  ISietNPAAI Tr 
SO  TO  A75 

continue 

LLPLt  NGIH<2) 

CALL  JNPACKIIOtlNtNUMAECtLL) 

00  223  A125  : It  NUXT-CC 

01  IM23  * lOFCOtTl  = 10  (K12Ei) 

CONTINUE 

lOFFGET  = lOFFS-lT  ♦ NUXSEC 

continue 

NAECIAA  : NAECIAA  ♦ NRECIN 
DO  -53  <150  : 1,  nj'<?EA 
e KISO)  = 01  (IGRIHA  ♦ <15.) 

continue 

HE  "EANS  ANO  THi  standard  OTVIATI'")''  of  TAE  TJO  -ATa  ATAAYS  AO- 
COMPUTED. 

SJ^KIA  I SU''<1R  0 SJNl<2A  o SUNTEO  = ' 

TNUNAEA  p NJAbEA 

OCNON  s FNUNBEA  • (FNJBBCA  - 1.1 
00  405  A152  = It  NUMBER 
SUNXlA  : SUNXIA  * A (R152> 

SUNXIB  = SUNXIB  ♦ B (1152) 

SUNX2A  P SUNX2A  * A (K152>  * A (K152> 

SUNX2B  : SUNX2B  * B IK152>  * B <K152I 
CONTINUE 

ANEAN  : SUMXIA  / FNUHBER 
BNEAN  = S'JNXIB  t FNUN3EA 

STOA  = IFNUMBER  » SUMX2A  - SUNXlA  • SUNXlA)  / DENOH 
STOA  = SORT  IA9S  (STOA>) 

STOB  = IFNUMBER  • SUHX23  - SUNXlB  • SUMXH)  / OENOM 
STOB  = SORT  (ABS  IST03)> 

THE  FORWARD  FFT  FOR  BOTH  DATA  ARRAYS  ARE  CONFUTED  USING  SUBROUTINES 
.REVBINt  CFFTRCt  AND  RTRAN2T. 

NUNEX*  s NUMEXF  * 1 
number  = 2 • • NUNExP 
FNJNBER  = NUMBER 
N = NUNEXP 
NM  = NUMEXP  - 1 
INV  = 1 
NDIR  = - 1 
SC  = .5 

CALL  REVBIN  (At  Bt  N) 

CALL  CFFTRC  (At  Bt  Mt  SCt  NDIR) 

CALL  RTRAN2T  (At  Bt  Mt  INV) 

THE  TRANSPOSE  Qt  THE  PRODUCT  OF  THE  FFTS  OF  THE  DATA  ARRAYS  IS 
CALCULATED. 

lU’LlN  : 1 ♦ NUMBER  / 2 
00  415  <160  r 1.  lUPLI" 

L160  = <16C  * lUPLlN 

real  = A (<160)  * A (L16b)  * B (<16G)  • B (L150> 

CONP  = B (K160)  • A (L160)  - A (<160  • P (L160> 

A (<160)  = real 
3 (<160)  = COMP 
CONTINUE 

THE  INVERSE  FFT  IS  TAKEN  OF  THE  ffT  PRODUCT  USING  SJ'ROUTINES 
REALTRAt  CFFTSt  AND  REORDER. 

INV  = - 1 
NOIR  = ) 

SC  = 1.  / FNUMRER 

CALL  REALTRA  (A.  Bt  MNt  NOIRt  INV)  » fc\l  f*  I 

C<LL  CFTTS  (At  Bt  MNt  SC.  NDIR)  . . f ill  |L|21  > 1 

CALL  REORDER  (A.  Bt  "N)  KM  fill  l\DL^  '■ 


- SUNXlA  • SUNXlA)  / DENOM 


- SUMXlB  • SUMXIB)  / OENOM 


B (<16G)  • 8 (LISO) 
A (<160  • P (L160> 


m 


KVMVA8«  tOPN 


imsn 
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15S 


160 


16S 


170 


175 


180 


185 


190 


195 


2oa 


C n:  LOCATION  AMO  i/ALJE  OF  Tit  OS  S-C9 « Af.  A TI  ON  TER*  IS 

C fOUNO,  ALONS  81H  Hi  '<IO=ni‘lT  OF  THl  31ST08TE3  0A1A  sroutNCE  HAT 

C LINfS  IT  UP  iJIH  THE  UNOTSTOSTEH  OATA  SEOJENCt. 

NUNhalp  = NUMBER  / ? 

NUMOTR  I NUMBER  / A 
'NJAHAL  = NUHHALF 
PNJMJTR  = NJMOTR 
PM£»CDR  : C. 

FCOR  = D. 

00  AS5  K170  = It  NU80TR 
FK170  s K170 

BLOUER  = B (K170  * NUNOTR)  / (FNUhOTR  * F<170  - 1.) 

01  (K170  * NUNQTR)  = IBLOWER  - A8EAN  • BMEANI  > STOA  / STOB 
RATIDI  = RFNUROTR  ♦ FK175  - 1.1  / FNU8HAL 

QUANl  s RATIOl  • 01  (K17U  * NUNSTR) 

IF  RQUANl  .L£.  FHAXCOR)  GO  TO  4:5 

FRAXCOR  3 QUANl 

FCOR  = 01  RK170  ♦ NUHSTR) 

LOCATIO  s K170  ♦ NUR3TR  - NUMBER  / 2 - 1 
425  BJPPER  = A RK170)  / IFNUNHAL  - FR17C  ♦ 1.) 

01  (K170  ♦ NUNMALFl  = (BUOPER  - ANEAN  • BREANJ  / STOA  / STOB 
RATI02  = RFNUNHAL  - FK17C  ♦ 1.)  / FNURHAL 
3UAN2  3 RATI02  * 01  (K17C  * NUNHALF) 

IF  RQUAN2  .LE.  PPAXCOR)  GO  TO  435 

FNAXCOR  3 0UAN2 

FCOR  3 01  <K170  ♦ NUHHALF) 

LOCATIO  3 K170  ♦ NUMHALF  - NUMBER  / 2 - 1 
435  CONTINUE 

C THE  LOCATION  ANO  VALUE  OF  THE  MAXIMUM  CROSS-CORRELATION  TERM  AND 

C THE  NEW  DISTORTED  TAPE  SEOUENCE  MIDPOINT  ARE  PRINTED  ALONG  ilTR  THE 

C INPUT  VARIABLES.  THE  NEV  DISTORTED  TAPE  SE3UENCE  MIDPOINT  IS  ALSO 

C PUNCHED  ON  A HOLLERITH  CARO  FOR  FUTURE  USE. 

PRINT  1502.MIDPNT1.'1I0PNT2.NUMEXP 
IFRMPRINT  .NE.  91  call  PRINTRStNUMPFB) 

PRINT  1522.  FMAXCOR.  FCOR.  LOCATIO 
1522  FORMAT  (2F12.3t  lOXt  15) 

MI 0PN3MI0PNT2 -LOCATIO 
PRINT  2,110. MIOPN 
2.:i0  FORMAT  «m) 

PUNCH  2010. MIOPN 
PRINT  2''00 
2')09  FORMAT  (3X.//) 

C THE  PROGRAM  IS  SENT  BACK  TO  READ  THE  NEXT  CONTROL  CARO. 

GO  TO  ns 
485  PRINT  1512 

1512  format  «.FND  of  PILE  READ  — JOB  TERMINATED  •) 

C PROGRAM  TERMINATION  SEJUENCF. 

4P5  PRINT  1514 

1514  format  (kHl.  • NORMAL  TERMINATION  •> 

475  CONTINUE 
END 


48 


BEST  mum  COPY 


s 


xo 


15 


20 


25 


I 


I 

I 

[ 


3C 


P<060»n  WRONIDP  (INPUTtOUTPUTtPUNCHX 
C PROGRAM  UR0MI3P  CALCULATES  THE  MIDPOINTS  Of  T^C  UOROS  OF  THE 
C DISTORTED  TAPE  KITH  RESPECT  TO  THE  UNOISTORTC)  TA«E. 

OINENSIDN  IPNISG) 

C THE  DRIFT  PER  SAMPLE  OF  THE  DISTORTED  TAPE  IS  CALCULATED  'ROM  •D|FF*t 
C THE  DIFFERENCE  IM  THE  NUMBER  OF  SAMPLES  BETWEEN  THE  UNDISTORTED  AND 

C distorted  tape*  and  •total**  the  TOTAL  NUMBER  OF  SAMPLES  IN  THE 

C DISTORTED  WORD  5R0UP  TAPE. 

READ  50*DIFF«TOTAL 
SO  FDRMAT(2F10.1> 

ORIFT=DIFF/TOTAL 
PRINT  35 
35  format  I1H1> 

C THE  50  UNDISTORTEO  PN  SEQUENCE  MIDPOINTS  ARE  READ  IN  FOR  THE 
C PARTICULAR  WORD  GROUP. 

DD  10  Isl,50 
READ  IS.MIPPN 
15  format  «I10» 

10  IPNfDsMIOPN 

C USING  THE  UNDISTORTEO  PN  SEQUENCE  MIDPOINTS*  THE  DISTANCE  BETWEEN 
C The  PN  SEQUENCE  AND  THE  WORD  OF  THE  UNOISTDRTED  TAPE.  AND  THE  DRIFT 
C PRESENT*  THE  WORD  MIDPOINTS  OF  THE  DISTORTED  TAPE  ARE  COMPUTED. 

00  2C  Jsl.EO 
REAO  25«NH0R0PN*N2S6 
25  format  II10*I5> 

IDRIFTfDRIFT.NHDROPN 

MIOWOROsIPNCJI.NWORDPNpIORIFT 

C THE  RESULTS  ARE  PRINTED  AND  PUNCHED  ON  HOLLERITH  CAROS  FOR  'UTURE  USE. 
PRINT  Sa.MIDU0R0.N256 
33  format  mo. 15) 

PUNCH  30.MIDU3RD.N256 
29  CONTINUE 
ENO 
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1^. 


£- — — ffgr*' 


1 


5 


IQ 


15 


20 


25 


SO 


35 


40 


♦ 5 


SO 


55 


' to 


05 


70 


C 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 


5 

30 


C 

C 

c 

c 

103 

105 


7 


C 


ID 

C 


25 

30 

15 


C 

c 

10 


*0 

*5 

S3 

55 

33 


MOOKAN  LOC  < II|>UT«OUT»UT«  TAPC 1*  Tt»e2s/t00l 

•ROCOtn  LRC  OOtS  TMC  7RANE  BV  7R«5£  L'C  PROCESSINS  Or  CACH  WORD  OF 
A PARTICULAR  MOlO  GROUP  AND  STORES  THE  IN'ORNATION  ON  NAGNETIC  TA^E. 
TNE  STORED  DATA  CONSISTS  OP  THE  ENERGT  TERR*  TRE  N0RMALI7E0 
AUTOCORRELATION  TERRS*  THE  PREDICTOR  COEPPICIERT  AUTOCORRELATION 
TERNS*  THE  PREDICTOR  COEPP IC lENTS*  THE  REFLECTION  COEFFICIENTS*  AND 

the  nininun  souareo  error. 

CONNON  A<25t)*ACOR(lS)*PCOEP(12*12>*6ANIlS>*ERROR 
OIRENSION  INCROI  « 10(490 1 *ITENPI9A‘!0),il(25tt 

INITIALIZE  VARIABLES.  *NUNBrR*  IS  EQUAL  TO  TRE  ANALYSIS  'RANE  SI7E* 
•NPC*  IS  EQUAL  TO  THE  NUHBER  OP  PREDICTOR  COEFFICIENTS  OESIREO* 

•ihoros*  is  the  nunber  of  sahples  in  a tape  recoro*  and  *>i*  is 

THE  UNIVERSAL  CONSTANT. 

NJRBERs25t 

NURHALFpNUNBER/2 

NPC»12 

IHORDSpROO 

TUsNJNBER 

»Is5.141592G535BB 

THE  HANNING  UINOOW  TO  RE  USED  ON  THE  SPEECH  SANPLES  IS  GENERATEO. 

DO  5 IsI.NUNHALF 

HH=(1.0/(I.0l»*TV>l*(0.S4*0.46*CaS(2.C*PI*I/TVI> 

W(12R*l>rW(12R-n2yH 

NAGNETIC  TAPE  IDENTIFICATION  VARIABLES  ARE  DE*1NC3. 

READ  3*NFILE*NTAPE*NTTPE 

FORHAT  (31101 

NRECORDsO 

NHORDsI 

PRINT  6 

FORNAT  (IHl) 

THE  NID“OINT  ANO  LENGTH*  IN  TERNS  OF  ANALYSIS  FRANE-,  OF  THE  VORO 
TO  BE  PROCESSED  IS  READ  INTO  THE  PRCGRAH.  IF  THE  VORO  HIDPOINT  IS 
EQUAL  TO  ZERO*  THE  ’ROGRAN  TERNINATFS.  INOEYES  ARE  C0N»UTE3  THAT 
POSITION  THE  tape  IN  ORDER  TO  GET  ^he  DATA  SAN»LES  FOR  TH£  JORO. 

READ  105*NIDPNT.N25S 
FORNAT  (110*15) 

IF  (N1D>NT  .EO.  0)  GO  TO  2>- 
NJLT=N256 

NLENGTH=NJNBER*N256 
PRINT  7,NI0PNT*NLENGTH*N256 
FORNAT  (2I1C.I5*/) 

INOEr=niopNT-NLENGTH/Z 

ISRIPsINDER/IWOROS-NRECORD 

I3RI5=INOE*-IWORDS*(INOE*/lyOROS) 

ALPHA  IS  defined*  VHICM  IS  USED  AS  AN  IDENTIFIER  ON  THE  NAGNETIC  TAPC . 
ENCODE  (19*h3«AL“HA>  NTYPE 
FORNAT  (17. SR) 

THE  TAPE  IS  POSITIONED  TO  THE  START  OF  THE  A0R3  TO  Dl  PROCESS'D. 

IF  (ISRIP  .LT.  1)  GD  TO  l' 

DD  ’5  KlSsl.TSRIP 

3JFFER  IN  (1*1)  (IN(1>.IN(m<')) 

IF  (IJNIT(1)>  15.29.25 
NPRRITYrNI'.NRECORD 


PRINT  39.NPARITY 

FORNAT  (*  PARITY  ERROR  IN  RECORO  **15) 

CONTINUE 

NRECORDsNRCCORD*ISKIP 

THE  SAMPLES  OP  THE  UORD  ARE  BUFFERED  INTO  THE  pROGRAN*  UNPACKED  BY 
SUBROUTINE  UNPACK*  ANO  STORED  IN  AN  ARRAY. 

NRCC=2*NLEN6TH/IMOROS 

IOFPSETsO 

DO  35  K35SI.NREC 

SUPPER  IN  (1*1)  (IN(1>*IN(B0I> 

IP  (UNITdl)  33*29*40 
KPARITY=K3S*NRCC0RD 
PRINT  4S.KPARITY 

fORNAT(*  parity  error  IN  RECORD  **15**  RUN  ABORTED  •) 

GO  TO  75 
LLpLENGTH(I) 

CALL  UNPACK  (ID.IN.LL) 

09  55  K55sl*IWOROS 

ITENP  (10PFSET4K55>xTD(K55) 

10PFSETsIOFFSET*IUOROS 
NRECOROsNRECORD*NRCC 
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TS  C ONE  Sr  ONE  THE  ANALYSIS  FAAHES  OF  THE  UOFO  AAE  JINDONED  BY  THE  HAMMING 

C WINDOW  AND  PMOCESSEO  BY  SUBMOUTINES  AUTOCOAt  *AEC3EF,  AND  GAMMA. 

DO  GC  KGOsltMULT 
NlMDEKsI0MI6*«MSfi-ll«N0M3EM 
00  AS  KASsItNUNBEA 

80  6S  AIK6S)sITeMP(KA5*NIN0E«)*U(K65) 

CALL  AUTOCOR  (NJM9EA«NPC> 

CALL  PRECOEF  INMCI 
CALL  GAMMA  (NPO 

C THE  L»C  INFORMATION  FOR  EACH  FRAME  IS  STORED  ON  MAGNfTIC  TA’E.  THE 

as  C first  six  VARIASLES  are  FDR  identification  purposes.  ACORIl)  IS  THE 

C ENERGY  TERM,  ACDRIJt,  J=?,13  ARE  THE  1?  NDRMAlIEED  AUTOCORRELATION 

C terms,  GAM(J),  J=1,1R  AFE  The  oREDICTOR  CDEFFICIENT  AUTOCORRELATION 

C TERMS,  »CaEFej,12»,  J=l,12  ARE  THE  OPEOICTOR  COEFFICIENTS,  9C0EF(U,JI, 
C J=l,12  ARE  THE  reflection  COEFFICIENTS,  ANO  ERROR  IS  THE  MINIMUM 

90  C SQUARED  ERROR. 

MRITE  I2,20'l)  ALPHA, NFILE, NT  APE, NMORO, MULT,  < AO,  (ACOR  ( J>  , J = I,  IT)  , < G 
1AM<  J>,J=l,n),(OC0EFIJ,12),Jil,12>,CPC0EFe  J,  J>,J=l,12>,FRRnR 
20  0 format  <A10,SIS,E16.13,12F12.10,13ElR.f',12ElR.R,12Fin.R,E12.A> 

S3  CONTINUE 

95  END  FILE  ’ 

backspace  2 

NHORDfNHORDmI 

NTY“E=NTypE*1 

C THE  PROGRAM  IS  SENT  3ACK  TO  READ  THE  NEXT  HORD  CONTROL  CARD. 

100  30  TO  1"0 

C PROGRAM  TERMINATION  SEQUENCE. 

20  PRINT  73 

70  format  (IMJ,«  normal  TERMINATION  ») 

T5  CONTINUE 

Its  END  FILE  2 

END 


1 
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15 


20 


25 
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PR06RAN  OISNE»$(INPUT«OUTPUT«T*PCl=/650tT*»C2=/CODtT«PE3=/lT5> 

C P106RON  OISNEAS  COMPUTES  TXE  RESIDUAL  OISTtNCi  MEASURES  D*«  0*  D*/E*« 

C AMO  D/E  AND  STORES  THEM  ON  MAGNETIC  TAPE  ALONG  WITH  THE  SUN  OF  THE 

C S3UARES  OF  THE  PREDICTOR  COEFFICIENTS*  THE  ENERGT  TERMS*  AND  THE 

C MINIMUM  S9UARE  ERROR  TERNS  OF  THE  UNOISTORTEO  AND  DISTORTED  TAPES 

C ON  A FRAME  NT  FRAME  BASIS. 

DIMENSION  AC0R(13>*GANI1S>*A(1?)*RC(12> 

DIMENSION  OACOR(13>*OGAM«13>*OA<12)*ORC(12> 

DIMENSION  E<*I*ASUNIB>*BC51*RI 

C INITIALI2C  VARIABLES.  *MUORn*  is  A COUNTING  VARIABLE. 

MW0RD=0 
PRINT  42 
42  FORMATtlHl) 

C THE  LENGTH*  IN  ANALTSIS  FRANCS*  OF  THE  WORD  TO  BE  PROCESSED  IS  READ 

C INTO  THE  PROGRAM.  IF  THE  VORO  LENGTH  IS  E3JAL  TO  ZERO*  THE  PROGRAM 

C JUMPS  OUT  OF  THE  PROCESSING  LOOP. 

50  REAO  10.N25G 
10  F0RNATII5> 

IF  IN256  .CQ.  0)  GO  TO  20 
HU0RD=NV0RD*1 
NMULTSN25S 
DO  38  1=1*8 
38  ASUMri)=a.O 

C frame  BT  frame*  THE  UNOISTORTEO  AND  DISTORTED  SPEECH  LPC  DATA  IS 

C REAO  INTD  THE  PROGRAM, 

00  15  K15=1*NMULT 

READ  ll*200>  ALphA*NFILE*NTAPE*NWORD.MULT*KSO*(ACOR( JI*J=l*13>t (GA 
IMI J)*J=1*13>«(AIJ>*J=1*12>*(RC( J) « J=1.12).ERR3R 
2G0  format  IAl?.5I5*ElG.10*12F12.1(l«13ElA.8*12ElA.8.12Fir>.8*E12.5l 

READ  <2*20'')  ALPHA.NFILE*NTAPE,NUORO,MULT*RSr,*<DACORIJ)*  J = l*13)  . ID 
1GAM( J)*J=1*13)*<DA( J>* J = 1*12)*(DRC( J> • J = 1 * 1 2 ) * DER R OR 
C THE  RESIDUAL  DISTANCE  MCASURFS  0.*  C,  D*/E**  AND  0/E  ARE  CALCULATED 

C AND  STORED  ON  MAGNETIC  TAPE  ALONG  WITH  SIA  IDENTIFICATION  VARIABLES. 

C THE  SUM  OF  THE  SQUARES  OF  THE  PREDICTOR  COEFFICIENTS.  THE  ENERGT 

C terms,  and  the  minimum  square  error  terms  of  the  UNOISTORTEO  AND 

C DISTORTED  TAPES  FOR  EACH  ANALYSIS  FR/iMr. 

SUM1=GAM<1) 

SUM2=DGAMI1) 

00  30  1=2.13 

SUM1=SUM1.GAM( l» .OACORI I) 

30  SJM2  = SUM2.0GAM(I  )*AC0R<I> 

E<1)=SJ»1 
E<2)=SU«2 
EI3)=EI1 )/DERROR 
E<A)=E<2)/ERRaR 
DO  32  1=1. A 

32  ASJMfI)=ASUM(l).C<I) 

ASUM<5)=ASUM<5)»AC0R(1) 

AS'J“»8)  = ASL'M<6>.0AC0R«  1 ) 

ASU"'<7)  = ASUM<7)*rRR'lR 
ASUH<  R)  = ASU'*(e)  ♦'lERROR 


/RITE  <3.3"!)  AL°HA.NFILE.NTAPE.NyORD.MULT.<fc3.E<l).:<2).E<5).E(‘) 
1.GAM(1).dg  AMU). ACnR)l).DACOR<l). ERROR. DERROR 
300  FORMAT<A10*3I5.1f.ElA.B) 

55  15  CONTINUE 

END  file  3 
BACKSPACE  3 

C THE  VORO  AVERAGE  OF  THE  RESIDUAL  DISTANCE  MEASURES.  THE  ENERGY 
C TERMS.  AND  THE  ERROR  TERMS  ARE  COMPUTED. 

GO  00  78  1=1.8 

78  8<MW0R0*I)=ASUM<  D/NMULT 
IF  <E0F<1))  20*25 
25  IF  <E0F<2))  20*A5 

C THE  PROGRAM  IS  SENT  BACK  TO  READ  THE  NEXT  JORO  CONTROL  CARO. 

G5  A5  GO  TO  53 

20  CONTINUE 

IFIMHORO  .EQ.  0)  MV0R0=1 
MM=NV0R0*1 

C THE  WORD  GROUP  AVERAGE  OF  THE  RESIDUAL  DISTANCE  MEASURES.  THE  ENERGY 
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TO 


7« 


RO 


85 
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C TtRISt  *ND  THt  £8908  TCR^S  *8E  CALCULATED. 

‘)0  90  1:1*8 
SU^sS.O 

DO  91  J=1«NW0R0 
91  SJ9sSU"*9CJ.I) 

9D  9l9U»HsSJ«>/9U39D 

C The  AOPD  average  and  the  word  grour  average  o*  the  residual  distance 

C HEASURES*  The  EVERG»  TERNS.  AND  THE  ERROR  TERRS  ARE  CRINTED. 

RRINT  T5 

T5  FORNATIJ'JH  COUNT  AR.A  A.RA.  El/E*  ET/E 

1 r:  erag  era6«) 

DO  SA  d:l,HUORr 
&R  RRINT  GS.J.cniJtl)*!:!*,  ) 

65  FORNATOIO.REIO.R) 

RRINT  51 
5S  FDR«AT«/> 

RRINT  6fi 

65  ROR«ATf»  VDRf)  S«OU»  AVERA5E  •) 

RRINT  £R«<9(NU.I  >«I:lt.E) 

AT  TDR^ATII  U*. 8^.10. *» 

C RRD3RAN  TERNlNATlDN  SEIUENC'  . 

RRINT  7 

70  FORMAT  <IH1,»  NORMAL  TERNINATICN  •! 

END  file  . 

END 


I 


i 
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PROGRAM  OISTKNdNPUTtOUTPUTfTAPCfs/U^tTAPf  3s/l  751 
PROGRAM  DISTKN  COMPUTES  THE  FACTOR  USED  IN  THE  •AI«  DISTAKCE 

MEASURE  AND  PRINTS  IT  FOR  EACH  WORD  AND  THE  UO> 0 GROUP  AVERAGE  FOR 
EACH  PAIR  OF  NOISE  AND  SIGNAL  FILES  READ  IN. 
dimension  EI2>tNU(5:>>FNAME(5> 

DIMENSION  F0AME(5I«0<2I «S(SOI 
DIMENSION  AC30)tR<3J>.OM(2>  tNCT(2> 

DIMENSION  Ql(51«5)«i2<51«5),33<51«5> 

DIMENSION  NQKSl  «5>,N02<51  .5  > .NCT.  <51  .5  ) 

INITALI7C  VARIA5LES.  .H'JORD*  IS  A COUNTING  VARIABLE. 

MWORD=0 

THE  LENGTH.  IN  ANALYSIS  FRAMES.  OF  THE  5t  WORDS  IN  THE  WORD  GROUP 
being  processed  are  read  INTD  THE  PROGRAM. 

Sf  READ  10.N256 
10  FORMATdS) 

IF<N256  .EO.  0)  GO  TO  2 ^ 

MWORD  = HHORO*l 
NMULT=N256 
NW<MUORO>=NMULT 
GO  TO  5 ! 

2t  CONTINUE 
HtfzMWORD*! 

A TOTAL  OF  .IFILl*  NOISE  AND  SIGNAL  DATA  FlLl'  ARE  READ  INTO  THE  PROGRAM. 
IFIL£=5 

READ  5.<FNAME<I) .1=1 .5> 

READ  5.(FDAMLd»  .I  = 1.5» 

5 FORMAT<5A7> 

EACH  OF  THE  PAIRS  OF  NOISE  AND  SIGNAL  DATA  FJIES  ARE  PROCESSED. 

DO  -5  R’Sil.IFILE 
LFN=LF0RM<FNAME<K25) } 

CALL  PFMATCH(5LTAPE J.LFF.r .O.O.i.l) 

LFN=LF0RM<FDAME<K25) ) 

CALL  PFMATCHCSLTAPED.LFN.r  . i.s.c.n 

EACH  OF  THE  WORDS  IN  THE  WORD  GROUP  FILES  ARl  prOCES'EU. 

DO  3 KI'.=l.MWOR!' 

N'>ULT=NW<K3.) 

SUM! =SUMi=U.O 
DO  -.5  I<15=1,NMULT 

FRAME  BY  FRAME,  THE  NOIif  AND  SIGNAL  LPC  DAtj  IS  read  INTO  THE  PROGRAM, 
read  <2.2''ni  ALOHA.NFILt  .NTAor  ,t,AORC.MULT,YfeU71,:2.ril>,0<2).Tl»T 
12.Y3.Y‘i.Y5.Yt 
DSC  F0RMAT<AlG,bI5»li1ElC.A> 

BEAD  <3.3ofc>  AL’HA.NFILE  ♦NTA-'f  .SWORE  .MULT.Kti  .yl.Y2,rd),E(2)<A:.n 
lAI, '-'■.dr:', ERROR. DERROR 
500  FORM'TCAIO, 515.1  ■EIR.m) 

ErtS  ‘R.  factor  IS  CD-PUTED  FOB  ' CCH  FRAME. 

A <.<  15  ) =«1 
S ( . 5 ) =7  ; 

B (■(:*  ):A 
SO“i=SU''l.l' J 

r.u’-'.  bSU'*;  ..!/’! 

5 CONTINUE 

HALE  the  AVEBa^C  tNf.-BGT  AND  thE  AVEtAGE  ‘K.  CTOR  O'"  THE  WORD 
BEING  PROCESSED  ARE  CALCULATED. 

RAWG=SUM1/NMULT 
KAVGI:RAVG/2.C 
ai<K3C.K2S>=SUM2 /NMULT 
Nal<R30.K25)=NMULT 

SUBROUTINE  HILOW  IS  CALLED  TO  COMPUTE  THF  A AGE  .R.  FACTOR  FOR  THE 

HIGH  AND  LOW  ENERGY  FRAMES  OF  THE  WORD  BEIN'  ROCESSEO. 
call  HIL0U<NHULT.RAVG2.A.R,S.0M.NCTI 
02<R3D.K25>=DHd  > 

Q3<K30.K25>=nM<2) 

Na2(X30.K25)=NCTIl> 

N03(K30.X25)=NCT(2> 

3(  CONTINUE 

The  average  factor  of  the  word  group  file  for  the  high,  low.  and 
total  Rnaltsis  frames  are  found. 

SUMlrSUM2=SUM3sO.C 
NCI rNCP=NC3=n 
DO  : 5 J=l.MWORO 
SU«USUM1«QI<J,R25>  , 

SUM, =SUM2.02( J,K2S> 

SUM* -SUM3.&3< J,R25) 
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NCl=NCl*N31(JtK!%) 

NC2=NC2*N02IJ«i<f5> 

NC5:NC3*NCi3CJtK2f>> 

23  CONTINUt 

01(Hm«K2  5I  :S(JNl/9UOf>D 
32«»’UtK25)sSUN2/"KORD 
33My,K2!.)=SUH3/''JO«ri 
tl&I  («Wt><2bl  :«iCl 
Na2C1M,K25»=*IC? 

\33l9y,K25)=fir.l 
?5  tONTlNUt 

C iHf  ’ iVESAbC  FACTORS  FOR  EACH  WORD  CF  EACH  WORD  GROUP  FILE 

C ARL  printed. 

PRINT  R 

RC  FORNATdHlt 
PRINT  R-.NFILE 
ur  FOR>AT<.  WORD  GROUP  NO, 

PRINT  7‘ 

73  FOR'  AT  ,M  VAR  I AT  LF 
1 ft 

CO  tR  Jrl.NJOaO 
PRINT  £=.J 

f>  FQR‘*iH.  wuRO  NUNRtR  ..I* 

PRINT  f 1 .(GI lU.I I flPlt' > 

F,.  ,<0.;  (U.I  > f I Pi  ) 
fci. (■).'<  J, I »,Isi»p| 
ht  » «N'.’  ( J.  1 >.  I PI  t‘  ) 

<-7.*N'...MJ.I)»IPl.'’  I 
Fp  < t )( J<  I > «lPl.N> 


.,I‘ ) 


GaiR 


19. 


9.F 


PRINT 
P.  IM 
-RI  NT 
.RIM 
I'k  I'.T 
PRINT 


1 1 •"  • 

35 

FOR' ATI/! 

NA 

CONTINUE 

G1 

FORMATdCH 

K 

.SEin.At 

62 

FORMATdOH 

KH 

•SEIP.A) 

63 

FORMATdJH 

«L 

tSriG.A) 

110 

86 

FORMATdOH 

NU^OCR 

,5110 

87 

FORMATduH 

NUN3£RH 

,5110) 

88 

FORMATd'iH 

NUHSCRL 

,5110 

115 


120 


12! 


PRINT 
PRINT 
the  3 
PRINT 


35 

35 

AVERAGE 

R5 


•<*  factors  for  each  word  group  file  are  printed. 


R5 


FORNATI*  WORD  GROUP  AVERAGES 
PRINT  GlttgilNU.IltlPl.*-) 
PRINT  G2t(02INW.I>«I=l*E> 
PRINT  G3,103CHW,n*I  = lR5> 
PRINT  86,(N31(NW,I>«lPl.5l 
PRINT  eT,(N32<NW«I).lPl.5l 
88*(Na3(NW«I>«Islf5> 
TERilNATION  SEQUENCE. 


• > 


7C 


PRINT 
PROGRAM 
PRINT  7 
FORMAT  OHl,. 
t NN 


NORMAL  TERMINATION  *> 


1 


PROGRAM  DISTSNnNPUT«OU1PUT»TAPC2i/l5%tTAP[3s/lTS> 

PROGRAM  OISTSN  COMPUTES  THE  LINEAR  DISTANCE  MCASURES  LMlHt  LMIL* 

LN2H<  and  LN2L  POR  EACH  UORD  AND  THE  WORD  GROUP  AVERAGE  FOR  EACH 
PAIR  OF  NOISE  AND  SIGNAL  FILES  READ  IN. 

DIMENSION  E(2)«NW(51>«FNAMCI^> 

DIMENSION  FOAMC(^>«D<2) .0(31,21 
DIMENSION  A(30,2>,R(30>.0M<4),NCT(2> 

DIMENSION  01(51.b),U2(SlfS>,93(Sl,5>,04(51,SI 

dimension  NOKSl  .b>.Na2(Ol,3> 

INITALIZE  VARIABLES.  *MUORO*  IS  A COUNTING  VARIABLE. 

MVOPOzO 

THE  LENGTH,  IN  ANALYSIS  FRAMES,  OF  THE  50  WORDS  IN  THE  UORD  GROUP 
BEING  PROCESSED  ARE  READ  INTO  THE  PROGRAM. 

51  READ  10,N256 
10  F0RMAT(I5> 

IF(N256  .EQ.  0)  GO  TO  2^ 

MUOROzMHORP.l 
NMULT=N2bG 
NH(MWOROI=NMULT 
GO  TO  5> 

2r  CONTINUE 
MWsMUORO*! 

A total  of  .IFILE.  NOISL  and  signal  DATA  FILES  ARE  READ  INTO  THE  PROGRAM. 
IFILEZ5 

READ  t,(FNAMC(I>,I=l,5l 
head  b,IFOAMEII> ,Iz1,5> 

5 F0RMAT(5A7> 

EACH  OF  THE  PAIRS  OF  NOISE  AND  riGNAL  DATA  FILES  APE  PROCESSED, 
no  .'5  RZbsl.b 
LFN:LF0RM«FNAME(K25I  > 

CALL  PFHATCH(5LTAPE3,LFN,:,  j,C,  ■,,!> 

LFN  = LF0RM<FDAHE<F25I  ) 

CALL  PFNATCH(5LTtPE2,LFN.u, 0.0,0, 11 

EACH  OF  THE  WORDS  IN  THi  WORD  GROUP  FILES  ARE  PROCESSED, 
no  3 K3C=l.HwORD 
N»ULTiN.(K3i  I 
•.U“  = t .0 

DO  15  K15=1,NHULT 

FRAME  BY  FRAME,  THE  NOISE  AND  SIGNAL  LPC  DATA  IS  READ  INTO  THE  PROGRAM. 
REAOIZ.ElO)  ALPHA,NFILC.NTAFE.NWORD.MULT,K60.Z1,72.D(1).C(2> ,Y1,Y2 
l,Y3,YR.Y5,Yi' 

2(0  FORMAT(Al->.tI5,l  ri  .A) 

ifAD  (3,300)  ALPmA.NFIL!  .•,TA>>F.NW0Rn.«ULT,K6l,*l,*2.E(l).E(2I.Ar.D 
1 AZ.RO.DRO.LHKCR.rERROR 
3iC  FOK''tT(A10.'l*.l  tlR.F) 

A(«:'..l)=ALOG(t(l)) 

A('(:b.2)zALf'G(E(.D>) 

(K :r.i) iALOG(D« , ) > 

‘•(Klb,2)zALDG(D(.>) 

R(F:b)=R; 

SUMrSUM.R,  • 

15  CCNTINUf 

•jALF  The  AYE'AGE  energy  op  T“E  word  '•FING  processed  is  COMPUTEn. 

RAVG2s(SUM/NMULT)/2.0 

SUBROUTINE  ENERGY  IS  CALLED  TO  COMPUTE  THE  DISTANCE  MCASURES  LMIH. 

LMIL.  LN2H.  AND  LH2L  FOR  THE  WORD  BEING  PROCESSED. 

CALL  CNCRGTINMULT.RAVGZ.A.B.R.DM.NCT) 
ai(K3I.R25)>DN(l > 

02(R30.K25)rDN(2> 

aS(K30.R25)s0M(3> 

Q4(R30.R25)zDN(«) 

N01(K3t.R25)sNCT(l» 

NO2(RS0.K25>sNCT(2> 

3r  CONTINUE 

the  r average  distance  measures  for  the  word  group  file  arc  found. 

SUH]sSUM2sSUM3zSUMRzO.S 

NClzNCZsC 

DO  2*  Jxl.NWORO 

SUMlsSUNI*ai(J.R25l 

SUM2=SUM2«a2(J.i(C5> 

SUM3zSUM3*03(U.'(:5t 

SUMRsSUN4«0R(J.<(2S> 

NClsNCl.NQK  J.K25) 


NC2zNC2.N92( J.R25t 


BESl  AVAIUBIE  COPY 


7b 


ao 


85 


9C 


95 


IGQ 


1(5 


110 


115 


120 


„.iii 


2* 


iS 


9P 

80 

75 


CONTINUE 

01INMtN2S>:$UNl/NU0R0 

a2CNUtK25>=SUN2/NU0KD 

03<MU«K25>sSUN5/HtfORO 

Q*l NUtK£5>  sSUNa/NUORO 

N01(HW«K05)=NC1 

N02<<(U*N25>=NC2 

CONTINUE 

THE  1VERA6C  DISTANCE  MEASURES  ROR  EACH  WORD  OR  EACH  NORO  GROUP  RILE 
ARE  PRINTED. 

PRINT  Rt 
RORHATIIHII 
PRINT  RC.NFILE 


OORO  GROUP  NO. 

variable 


69 


RORNATt* 

PRINT  75 
FORNATCbOH 
1 AR» 

DO  6A  Jsl.HRORO 
PRINT  69, J 

FORMAT**  VORO  NUHnCR 


iI5) 


56.9 


19. 


9.6 


FN 


*.1*) 


PRINT  61 

J,l >,I: 

1*5) 

PRINT  62 

f <r.^( 

J*I >*IP 

1,5) 

PRINT  63 

Jtll.lR 

1«5> 

PRIM  66 

• (Gi4( 

J,I >.IA 

!»5> 

PAINT  86 

• <NU1 

*J, I >,I 

Pl,5» 

PRINT  67 

f iNor 

*j.t>*: 

= 1.5  > 

PRINT  35 

35 

FORMAT*/ 

) 

69 

CONTINUE 

61 

FCR«iT*l 

L»1H 

. 

« 1 

FORMAT*! 

'H 

l»;h 

.5C i ' . 

b5 

FOR  ■•AT  *1 

LMIL 

.5f .£. 

4 > 

FORMATIl 

OH 

LM2L 

•5Eie. 

4> 

66 

RORMATEl 

OH  NUMBERH 

•51101 

67 

FORMATdOH  NUHBERL 

*5110) 

95 


70 


PRINT  35 
PRINT  35 

THE  AVERAGE  DISTANCE  NEASURCS  FOR  EACH  WORD  GROUP  RILE  ARE  PRINTED. 
PRINT  95 

FORNATt*  HORD  GROUP  AVERAGES  •> 

PRINT  61*iallNW*I>«Isl*5) 

PRINT  6P*(Q2INH«I)*I:l«5> 

PRINT  63*(tl3INU«II«lPl«5l 
PRINT  66»tQ9{HV«I>*lzl«5) 

PRINT  S6«CNQllHW*II*Izl,5l 
PRINT  a7*(NQ2<NU*I>*lPl,5) 

PROGRAN  TERMINATION  SEQUENCE. 

PRINT  Tr 

FORMAT  (IHlt*  NORMAL  TERMINATION  *1 

END 


BEST  AVAILABLE  COPY 


57 


vrx..r-:. 


r 


SU4<0UTINt 
SUSWOUTIDE  PlIIKn 


IN*  NliMttP) 

p^iMTS  *tL  or 


10 


IS 


Zv 


25 


ns 


120 


CALCULATED  BT  F^TCOAA  FQA  A PARTtCULAA 
CORSON  /FFTl  /A  r*20nit  B (A200t«  3UN1 
lOl 

CONNON  /FFT2  FICOR  (210S> 

NUNHALF  s nurse a / 2 
NJNSTR  z NUNSER  / A 
«U»LIN  » 1 ♦ NJSOTR  / 1C 
LMIT  X NUNQTR  - 1 
C (LtNtT  » 1)  X 0. 

DO  115  K2S  X It  LIPIT 
C CK25>  X 0 

C IK25  * NUNHAL^  • NUNQTRI  x q 
CONTINUE 

DO  120  KSO  X 1,  KUPLIN 
LCNTR  x - 10  • IKUPLtS  - KSOI  - 1 
NCNT*  X LCNTA  - 9 
LL  X NUSHALF  ♦ 1 ♦ LCNTR 
PRINT  ISOOt  LCNTR*  C ILLit  C (LL  - Dt  C 
ILL  • Alt  C (LL  • 51t  C (LL  • 6lt  C (LL  • 
2>t  ACNTR 
CONTINUE 

PRINT  1502t  C (NURNALF  ♦ 11 
30  125  RAC  X It  FJPLIR 


the  cross-cdrrelation  t<:rns 


SET  DF  NIDpOINTS. 

(A2(3>t  DUN2  (A2J0>t  C (SCO 


(LL  - 2>t 
Tit  : (LL 


(LL  - *1,  C ( 
Bit  : (LL  - 9 


30 

LCNTR  X 10  • (KA3  - 11  ♦ 1 

NCNTR  X LCNTR  ♦ 9 

LL  X NUHHALF  . 1 ♦ LCNTR 

LJP  X LL  ♦ 9 

PRINT  15C0.  LCNTR,  (C  (I>,  I x LL,  L'J’l.  NCNTR 

1 

125 

CONTINUE 

1500 

RETURN 

FOR«AT  (15,  2Kt  1CF12.5.  2»t  151 

1512 

FORNAT  (.  t .,  12*.  PI?. 31 

.y 

55 

END 

1 


s 


10 


»9 


20 


SUBROUTINE  UNPACK  (lOtINtLL) 

C SUBROUTINE  UNPACK  TAKES  A SO  BIT  MORO  AND  SEPERATES  IT  INTO  FIVE 
C 12  SIT  SANPLESt  THIS  IS  DONE  FOR  THE  EACH  SC  HORO  RECORD. 
OINENSION  IO(A00ltIN(80l 
LSx>3 

DO  SB  IxltLL 

LT«IN(II 

LSxLS*lB 

DO  SO  L«ltS 

LTIxlT  .and.  S77TB 

ltjxlt  .and.  AOOOB 

IF(LTJ  .NE.  01  LTIx  .not.  LTI 

IDtLSlxLTI 

C PSNIFT(LTt'12l*  SHIFTS  THE  BITS  OF  VARIABLE  xLT.  TO  THE  RI6HT  12 
C ' ' •(.sets.  THIS  IS  PERFORHEO  IN  ORDER  TO  OBTAIN  THE  FIVE  12  BIT 
C SSNPLES  FROH  each  S'*  BIT  HORD. 

Lr«SHlFT(LTt-12> 

LSxLS'I 
SB  CONTINUE 
SO  CONTINUE 
RETURN 
END 


..  BEST  AVAIUBIE  COPY 


1 


s 


10 


IS 


20 


25 


30 


35 


«0 


*5 


SUBROUTINE  REVBIN  («t  B»  HNI 
C CALL  REVBINfR*B*H» 

C REVERSIBLE  PERMUTATION  OP  ARRAYS  A AVD  B 

C PRON  N3RNAL  SEQUENCE  TO  REVERSE  BIVART  SEQUENCE* 

C OR  VICE  VERSA* 

C SEQUENCE  LENBTH  TS  N * 2**N 

C VRtTTEV  BY  L*  OAVIO  LEVIS  AND  NARIE  VEST*  ESSA. 

C NOOIPIEO  FROM*  OR  INSPIRED  BY  THE  AL30L  PROCEDURE 

C REVERSEBINART*  BY  R.  C.  SINSLETON*  SRI. 

DINENSIDN  A fltSBAt*  B IlSSRA) 

CONNON  /FFTCC  /M*  UC  (ISI*  ST  (151 

>1  X AN 

CALL  PPTC 

IF  (N  .LE.  II  RETURN 
N X JC  (N  * 11 
HP  X N ♦ 1 
R X 1 
I X 2 
J X N - 1 
IfO  LC  X N 
105  R X K ♦ JC  (LC) 

JC  (LC)  X . JC  (LC) 

IF  (JC  (LC)  .LT.  0>  SO  TO  110 
IF  (LC  .EQ.  2)  RETURN 
LC  X LC  • 1 
SO  TO  1C5 

110  IF  (K  .LC.  1 .OR.  J .LT.  K>  GO  TO  115 
T X R (I> 

A (I>  X ( (K) 

A (R)  X T 
T X 3 (1) 

9 (I)  X R (K> 

J (K)  X T 

IF  (J  .E9.  R>  63  TO  115 
RR  X N»  - R 
T : A (RR) 

A (KK>  X A (J) 

A (J)  X T 
T X 3 (RRI 


R (RR) 

S ' 

3 ( J) 

= T 

I X T 

♦ 1 

J X J 

» Q 

63  TO 

I'O 

END 

BEST  AVAIUBIE  COPY 


SUSROUTINE  CFFTRC  <«•  Rt  RRt  SC«LCi  NEXF> 

DISCRETE  CONPLEX  FAST  FDJRIER  TRANSFORM. 

CALL  C^FTRCIA«B*n»SCtNX) 

INF'JT  A(J>  * I*P(JI  IN  REVERSE  BINART  SEQUENCE. 

OUTPUT  ACK>  * I*SIK>  IN  NORMAL  SEQUENCE. 

SEQUENCE  LENSTH  IS  H : Q**N 
SC  IS  REAL  SCALIN6  MULTIPLIER. 

NX  IS  THE  SI6N  OF  THE  EXPONENT  IN  THE  TRANSFORM  DEFINITION 

INNER  LOOP  SINES  AND  COSINES  COMPUTED 

RECURSIVELT  3V  SINGLETONS  ?NO  OI^'ERENCE  ALGORITHM. 

initialized  from  a data  table. 

URITTEN  BT  L.  DAVID  LEVIS  AND  MARIE  WEST.  ESSA. 

MODIFIED  from.  OR  INSPIRED  AT  THE  ALSDL  PROCEDURE 
REVERSEFOURIERC.  BV  R.  C.  SINGLETON.  SRI. 

DIMENSION  A IlSBBAI.  B 116384) 

COMMON  /FFTCC  /M.  JD  flS).  S (15) 

M s MM 

CALL  FFTC 

N s JD  (M  * 1) 

K = N / 4 

NO  s K 
NN  s N - I 
JSPAN  : 1 
SC  s SCALE 

IF  (ABS  (SC  - 1.)  .LT.  l.E-ID)  GO  TO  IDS 
DO  iro  JC  .=  1.  N 
A (JC)  s SC  • A (JC) 

B (JC)  = SC  • B (JC) 

108  CONTINUE 
105  IF  (M  .EQ.  S)  RETURN 
DO  110  MM  = 1.  N.  2 
MS  s MM  ♦ 1 
RE  s A (MM)  - A (MS) 

A (KM)  = A (MM)  * A (MS) 

A (MS)  s RE 

=1M  = B (MM)  - B (MS) 

R (KK)  : 3 (MM)  * B IMS) 

9 (MS)  = FIM 
118  CONTINUE 

IF  (M  .tO.  1)  RETURN 
EXPS  = 1S16N  II.  NEXP) 

N»  s 1 

DO  125  J9  = 2.  M 
SD  = - S (JB  - 1) 

CO  = 2.  • S (JB)  • S (JB) 

R = - 2.  • CD 
CN  = 1. 

CM  = 0. 

SN  r 0. 


SM  s ♦ r»“S 

JSP4NH  I JSPAN  _ 

JSPAN  s JSPAN  « JSPAN 
115  KS  s KK  * JSPAN 

RE  3 CN  • A (KS)  - SN  • B (KS) 

FIM  s SN  • A (KS)  * CN  * 3 (KS) 

A (KS)  3 A (KK)  •>  RE 

A (KK)  3 A (KK)  * RE 

B IKS)  3 B (KK)  - FIM 

3 (KK)  K 3 (KK)  * FIM 

KK  3 KK  • JSPANH 

KS  3 KS  * JSPaInH 

fin  3 SM  • a (KS)  « CM  * B (KS) 

RE  3 CM  • A (KS)  - SM  * B (KS) 

) <;.♦  R * .fi"*  - "E 

A,IKK»  3 A ;(KK)  ♦-  RE 
3 (KS)  3 B (Kk)  • FIM 

B (KK)  3 B (KK)  ♦ FIM 

KK  3 KS  * JSPANH 
IF  (KK  .LT.  N)  SO  TO  115 
KK  3 kk  - NN 

JJ  3 JJ  * K 

IF  (JJ  .SE.  NB)  60  TO  120 


75 


SC 


35 


C3  = R • CN  ♦ CO 
CN  r CO  ♦ CN 
’IH  = CN  • EX  = S 
I R • CM  ♦ S3 
CM  = SO  ♦ CM 
SN  = - CM  . £«os 
GO  TO  115 
120  M : K / 2 
125  COMTIMUi 
RETURM 
END 


1 


5 


10 


15 


20 


25 


33 


35 


RC 


R5 


5u 


C 

c 

c 

c 

c 

c 

c 

c 

c 


130 

105 


110 


115 

123 


SUBROUTINE  RTRAN2T  U.  MM<  INVl 
CALL  RTRAN2TIAfRtMf INVl 

IFIINV.GT.P)  UNSCRAMBLE  THE  TRANSFORMS  OF  TWO  REAL  SEOUENCES. 
IFIINV.LT.OI  scramble  the  transforms  3F  two  real  SE3UENCES. 
INPUT  ANO  OUTPUT  ARE  IN  NORMAL  SE3JCNCC< 

SEE  WRITEUP  FOP  DETAILS. 

SEQUENCE  LENGTH  IS  N = 2*«M 

WRITTEN  BY  L.  DAVIO  LEWIS  ANO  MARI-  WEST.  ESSA. 

MODIFIEO  from,  or  INSPIRED  PY  THE  ALGOL  PROCEDURE 
REALTRAN,  BY  R.  C.  SINGLETON,  SRT. 

DIMENSION  A (1&3R5>,  3 tl&3HGI 
COMMON  /FFTCC  /M,  JC  «1E),  ST  <151 


M s mm 

call  fftc 

N = JC  (M  ♦ 11 
NH  s N / 2 

IF  UNV  .LT.  01  GO  TO  120 
IF  IM  iLT.  21  GO  TO  105 
K = N 

DO  i:0  J = 2,  N-l 
A (K  * 1)  = R (J>  M P CKI 

R IK  ♦ 11  = A IMI  - A (J) 

A (Jl  : A CJI  * A (K) 

3 <J>  = B (Jl  - B (K) 

K s < - 1 

A 111  3 2.  • A (II 
A <NH  *21  = 2.  • R <11 
3 (1)  s 1 (NH  ♦ 21  I 0. 

IF  <M  .E3.  0)  RETURN 
A (NH  * 1)  = 2.  • A (NH  ♦ 11 
A (N  ♦ 2)  = 2.  • R (NH  ♦ 3) 

3 (NH  ♦1>=3(N«2>=0. 
IF  (M  .EQ.  2)  RETURN 
< z N ♦ ' 

J z NH  ♦ 3 
NS  = NH  / 2 - 1 
DO  115  L = 1,  NS 
T z A (Jl 
A (Jl  z A (Kl 
A (<1  z T 
T z R (Jl 
3 (Jl  z B (Ml 
P.  (K)  = T 
J z J « 1 
R z K - 1 
RETURN 

3 (11  = A (NH  ♦ 21 

IF  (M  .LG,  01  RETURN 


BEST  AVAIWBIE  COPY 


^ (NH 
1?  (w 
< s NH  ♦ 2 
DO  121  J = 

8 (Kl  s A 


1)  s A (N  ♦ 2) 
iQ.  1>  «»£TjqN 


2*  NH 
( J) 


^ (< 


I ) 


61 


95 


6Q 


a <Ki  s * (K  * 1 ) - e ( j) 

« (j>  s « (j)  • a (K  * 1) 

3 IJ>  s B IJ>  » A IB  » II 

125  a = K ♦ 1 

IF  «a  .EO.  2)  RETURN 
R r N 

J = NX  ♦ 2 
GO  TO  lie 
END 


i SUBROUTtNF  FFTC 

C CONNON  SUBROUTINE  FOR  FFT  SUBROUTINES. 

C JC  IS  »OWERS  OF  TWO  ARRAY..  JC (N I : 2 •* (M-1 I 

C ST  IS  SINE  ARRAY..  ST ( n)sS INI  PI /f 2 • •1) I 

5 C B IS  TESTED  FDR  PROPER  INPUT  RAN3E.  I.LE.H.LE.IA. 

CONHON  /FFTCC  /N.  JC  1151.  ST  115) 

DATA  IJC  : 1.  2.  A.  B,  16.  52.  6A.  123.  25G.  512.  1G2A.  20*3.  AGSG 

1.  B192.  1658*1 

OAT*  (ST  s l.OOOOOOOOOOCE^OOO.  7. 0 71C 6781 1 B7E *00 1 . 5.82665*525651- 
10  1001.  1.9S0NO522016E-001.  9.80171*052961-002.  *.90&767*327*C-002.  2 

2. *5*12235229E-092.  1 .22715382857E-0 92.  6.15583*6*9151-005.  5.96795 
3676297E-005.  1 .55598018629E-C05.  7.6699''31  3 7*5 E-00* . 5.35*95187571 
*E-00*.  1.917A7597311E-00*.  9.537579909691-005) 

IF  (>*  .LT.  0 .OR.  N .GT.  1*1  PRINT  1500.  7 
15  RETURN 

1500  FORNAT  KCILLESAL  VALUE  FOR  N«  N =*.  1*1 
ENO 


1 


SUEKOUTINC  REALTR*  It,  B,  NC « IRV) 

CALL  REALTRASIA«B«R*NC*I*IV> 

IF(INV.GT.3)  UNSCRAMBLE  THE  TRANSFORM  OF  A REAL  SIOUENCE. 
IF(INV.LT.D>  SCRAMBLE  TMC  TRANSFORM  OF  A R£*L  SEOUEMCE. 

INRUT  ANO  OUTPUT  ARE  IN  NORMAL  SEQJENCEt 
SEE  URITEUR  RQR  DETAILS. 

SEQUENCE  length  IS  N s 2..H 

NE  MUST  AGREE  WITH  SIGN  OF  EXPONENT  IN  TRANSFORM  DEFINITION. 

INNER  LOOP  SINES  ANO  COSINES  CO*»UIEO 

RECURSINELT  ST  SINGLETONS  2NO  OIF^ERiNCE  ALGORITHM, 

initialiteo  from  a data  table. 

ARITTEN  ST  L.  DAVID  LEUIS  AND  MARIE  VEST,  ESSA. 

MOOIFIEO  FROM,  OR  INSPIRED  ST  THE  ALGOL  PROCEDURE 
REALTRAN,  ST  R.  C.  SINGLETON,  SRI. 

DIMENSION  A I1AS85>,  3 (16S8SI 
COMMON  /FFTCC  /M,  JC  fl5),  ST  <15) 

M = MM 
CALL  FFTC 

N = R = JC  ♦ I> 

NH  = N / 2 

NX  = NH  ♦ 1 

CN  = ISIGN  (1,  INV) 

SN  : ISIGN  (I,  NE) 

A «NK>  : 2.  . A INK) 

3 INK)  = 2.  . CN  . SN  • R CNK) 

IP  <CN  .GE.  0.)  60  TO  iro 
PIM  = A (1)  - A (N  . 1) 

All>  = All)»Ar)«l> 

3 <n  = FIM 
GO  TO  1-5 

A IN  ♦ 1)  = 2.  • CA  ID  - R ID) 

A ID  = 2.  . lA  ID  . R ID) 

3 ID  = a IN  ♦ D = n, 

IF  IM  .ED.  j)  RETURN 


SD 

= CN  • SN  * 

ST 

= 2.  . ST  (H 

* 1) 

= - R . R 

CO 

= - .5  . CN 

• 3 

SN 

: o. 

'^•0 

MO  J = 2, 

N-t 

CO 

i R » CN  ♦ 

CO 

CN 

= CO  ♦ CN 

SO 

= R . SN  ♦ 

SO 

SN 

r SO  ♦ SN 

ftA 

= A t.t)  * A 

( Kl 

ft9 

= A IJ)  - A 

(<) 

3A 

: a IJ)  ♦ P 

<•() 

33  : 3 IJ)  - 3 IK) 

RE  = CN  . 3&  ♦ SN  . A3 


50  PI“  - SN  . BA  - CN  * AR 


SUBROUTINE  CEFTS  C*t  9i  HN«  SCRLE*  NElR) 

DISCRETE  C3RPLER  FAST  FOURIER  TRANSFTRR. 

CALL  CFFTSIA«B<NtSCtNX> 

INPUT  A(UI  * I*B(UI  IN  NORNAL  SESJENCE. 

OUTPUT  AIK)  » I*8IN)  IN  REVERSE  BIRARV  SEOJENCE. 

SEQUENCE  LENSTH  IS  N s 2**N 
SC  IS  REAL  SCALING  NULTIPLIER. 

NX  IS  THE  SISN  OF  THE  EXPONENT  IN  THE  TRANSFORM  OEpINITION. 
INNER  LOOP  SINES  AND  COSINES  COMPUTEO 
RECURSIVELY  ST  SINGLETONS  2ND  DIFFERENCE  ALGORITHM. 
INITIALIEEO  FROM  A DATA  TABLE. 

VRITTEN  RY  L.  DAVID  LEVIS  AND  MARIE  VEST.  ESSA. 

MODIFIED  from,  or  INSPIRED  RT  THE  AL30L  PROCEDURE 
fastfojriers.  By  r.  singleton,  sri. 

DIMENSION  A IlGSRR).  3 116384) 

COMMON  VFFTCC  /M.  JC  115).  SNT  115) 

MA  = M = MM 
CALL  FFTC 
N = JC  CM  ♦ 1) 

NH  = N / 2 
N1  r NH  / 2 

SC  = scale 

IF  <M  .El.  0)  30  TO  120 
IF  |«  .E2.  1)  30  TO  110 
EX9S  = ISIGN  (1.  NETP) 

NN  = N - 1 

K I ; 

30  110  JA  = 2.  H 
CE  = SNT  IMA) 

MA  = MA  - 1 

CD  r 2.  • Ce  • ZC 

SD  : - SNT  (MA) 

R = - 2.  • CD 
CN  = 1. 

CM  = 0. 

SN  = 1. 

JJ  = "1 

FK  = 1 

SM  = ♦ E»PS 

JS’AH  = NH 

NH  = JSOAN  / 2 

KS  = HH  ♦ JSPAN 

RE  = A (<K)  - A (KS) 

A (KK)  = A (KK)  ♦ A (KS) 

=IM  : B (KK)  - H (KS) 

3 (KK)  = 3 (KK)  ♦ 3 (<S> 

A (KS)  I CN  • RE  - SN  • FIM 

3 (KS)  : SN  * RE  ♦ CN  • FIM 

KK  I KK  ♦ NH 
KS  = KS  ♦ Nh 
RE  = A (KK)  - A (KS) 

A (KK)  I A (KK)  » A (KO) 

PIM  r R (KK)  - 3 (KS) 

B (KK)  A B (KK)  * M (KS) 

A (KS)  = CM  . RE  - SM  « FIM 

8 (KS)  = SM  . RE  ♦ CM  . FIM 

KK  = KS  * NH 
IF  (KK  .LT.  N)  30  TD  ICO 
KK  = KK  - NN 
JJ  = JJ  • K 


IF  (JJ  .GE.  NQ)  GO  TO  1S5 

CD  = R K CN  * C3 

CN  = CO  ♦ CN 

SD  = R • CM  ♦ S3 

CM  * (i»t  ♦ SD 

SN  z ^ CM  • EKPS 

SM  z CN  • EKPS 

SO  TO  110 


K = K * K 
CONTINUE 

DO  120  KK  = 1.  N.  2 

KS  = KK  . 1 

RE  = A (KK)  - A (KS) 
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* «KS)  s RF 

F11  = B «KK>  - H IKS) 

3 IKK)  : (B  IKK)  * 3 IKS)) 

B IKS)  : FIN 
CONTINUE 

IF  lABS  ISC  - 1.)  .LT.  i.E-lO)  GO  TO  135 
00  130  JB  s If  N • 

* IJB)  : SC  f * US) 

B IJB)  = SC  • 3 IJB) 

CONTINUE 

return 

END 


SUBROUTINE  PRECOEF  INPC) 

SUBROUTINE  PRECOEF  CONPUTES  THE  12  PREDICTOR  C9EFF IC lENTSt  THE  12 
REFLECTION  COEFFICIENTSf  *ND  THE  NININUH  S3JRRED  ERROR  USING  THE 
LEVINSON  algorithm  AND  THE  AUTOCORRELATION  TERNS  FROM  SUBROUTINE 
AUTOCOR.  SUBROUTINE  STABLE  IS  USEO  TO  INSURE  A STABLE  FILTER  IS 
generated.  the  predictor  COEFFICIENTS  ARE  STORED  IN  »C0EFUfl2)f 
J=lfl2t  the  reflection  coefficients  are  stored  in  PCOEFFIJtJIf 
J=lfl2f  AND  THE  NININUH  SDUAREO  ERROR  IS  THE  VARIABLE  0. 

C0NN9N  AI256)f ACORI13)fPCOEFI12f 12)fGANI 13),D 
C=-AC0RI2) 

D=1.0“C«C 
PCOEFIlf 1)=-C 
□0  10  J=2fNPC 
J1=J-1 
SUNin.o 
DO  15  IslfJl 

SUNsSUNfACORIJ*l-I )fPCO£F|l,Jl ) 

GsACORI J*1)-SUN 
C=-G/3 

IF  IA3SIC)  .GT.  0.99)  CALL  STABLE  IJ.C) 

D=D»I1.0-C‘C) 

DO  20  I=lfJl 

9C0EF|I,J)=PC0EFlI,Jl)fCfPC0EF|J-IfJl) 

’CDEFUfJ)=-C 

RETURN 

END 
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85 


120 

125 


130 

155 


10 


15 


20 


25 


15 


20 

10 


1 


c 

c 

c 

s c 

c 
c 
c 

la 


IS 


suaaouTiNe  acoaocii  la*  St  •«(> 

CiLL  SEOSOERia.Bt'O 

REVEKSISLC  BEBBUTaTION  OF  KEAL  SEO'JE>tCE 
FROB  FIRST-LRST'NORNRL  SEQUENCE 
TO  OOO'EVCN-RCVERSE  BINARY  SEQUENCE. 

SEQUENCE  LENQTH  IS  N s 2**N. 

WRITTEN  BY  L.  DAVID  LEVIS  AND  MARIE  VEST.  ESSA. 
REORDER*  9Y  R.  C.  SINGLETON*  SRI. 

DINCNSION  A «1GSBA>*  9 tlGSBR)*  LST  (IS) 

COMMON  /FFTCC  /M*  LC  l»5)*  ST  (15) 

M X MM 
CALL  FFTC 

IF  (N  .ED.  0)  RETURN 
JA  X M ♦ 1 
JB  X M.-  I 

1 X RB  X !) 

NU  X LC  (JA)  - 1 
DO  ICO  KA  X 1,  KU.  2 
T X A (KA  * 1) 

A (KA  * 1)  X B (KA) 

3 (KA)  X T 
CONTINUE 

IF  (N  .EQ.  1)  RETURN 
LIN  X M / 2 * 1 
KS  X LC  (JB  * 1)  * KB 
KU  X KS 

JJ  X LC  (JA  - J3) 

KK  = KB  * JJ 
K X KK  * JJ 
r X A (KK  * 1) 

A (KK  t ] ) X A (KS  * 1) 

A (KS  ♦ 1)  X T 
T X B (KK  * 1) 

X (KK  * 1)  X B (KS  ♦ 1) 

3 (KS  ♦ 1)  X T 
»K  X KK  ♦ 1 
KS  X KS  ♦ 1 

IF  (KK  .LT.  K)  iO  TO  11? 

KK  X KK  * JJ 
KS  X KS  ♦ JJ 

IF  (KK  .LT.  KU)  SO  TO  IIQ 
IF  (J9  .LE.  LIN)  SO  TO  :2C 
J3  X JB  - 1 
! X I ♦ 1 
LST  (1)  X J3 
50  TO  IPS 

1'  (1  .LE.  j)  RETJRN 
JB  X LST  (1) 

1 = 1-1 
K3  X KS 
SO  TO  1J5 
END 
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SUBROUTINE  RUTOCOR  (NUNRCRtNPCI 

SUBROUTINE  AUTOCOR  CONFUTES  THE  ENERST  TERN  AND  THE  IZ  NORHALIEEn 
AUTOCORRELATION  TERNS  FOR  EACH  2SB  SARBLE  Ml  NOOK.  THESE  VALUES 
ARE  STORED  IN  ACOR<JI«  RESPECTFULLY. 

CONNON  A<25A>*AC0RllS)tPC0EFI12*S2>t6ANCI5>*CRR0N 

NsNPC*l 

SUN«B.0 

00  10  Isl.NUNBER 
13  SUNsSUN.ACII 

SNEANsSUH/NUNBER 
00  15  Isl.NUNBER 
15  AfItsAIIl-SNEAN 
SUNsO.O 

30  20  Irl.NUNSER 
20  SUNsSUN*ACII*Altl 
ACORdlsSUN 
DO  25  JsZ.N 
S'JNcO.P 

RUNsNUNBER-J*! 

00  30  1*1 .NUN 
SO  SUNsSUN.A(1>*AII*J-1> 

25  ACOR(J>sSUN/ACOR(l> 

RETURN 

ERO 


SUBROUTINE  STABLE  fUtCi 

SUBROUTINE  STABLE  FORCES  THE  ABSOLUTE  VALUE  0*  EACH  REFLECTION 
COEFFICIENT  TO  BE  LESS  THAN  OR  E9UAL  TO  B.BB*  THFREBT  INSURINS  A 
STABLE  FILTER. 

CCa-C 

PRINT  lO.J.CC 

ID  FORNAT  I.  REFLECTION  COEFFICIENT  •*I2t*  IS  UNSTABLE.  « •t*lD.S> 
IF  1C  .BT.  O.BBl  CaO.BB 
IF  CC  .LT.  -O.BBt  Cf-B.VB 
RETURN 
END 


SUBROUTINE  BANNA  IN*C> 

SUBROUTINE  BANNA  CONFUTES  THE  PREDICTOR  COEFFICIENT  AUTOCORRELATION 
TERNS  FRON  THE  PREDICTOR  COEFFICIENTS. 

CONNON  AI25BltAC0R(lS>«PC0CFI12f 12) *SAN( 131 .ERROR 
OINENSION  PCfl3) 

NaNPCl 
PCIllPl.O 
DO  10  I«1«NPC 
10  PCII«l>>-PC0CFfl«12) 
suNso.a 
DO  15  Isl.N 
15  SUNaSUN«PC(l).P:Cl> 

SANIllsSUN 
DO  20  Isl.NPC 
SUNsO.O 

„ BEST  AVAILABLE  CO 

20  SANri«l)«2.0*SUN 
RETURN 

END 


1 


c 

c 


f 


1C 


19 


2ft 


SUBROUTINE  HILOU(NMULT«R «V62tA«R •S*ON«NCT) 

SUBROUTINE  HILOM  COMPUTES  THE  AVCRROr  (K*  FACTC»  FOR  THE  HI6H  AND 
LON  ENCRST  FRAMES  OP  THE  UORO  BEING  FROCESSED* 
niNENSION  AISOltRISMtSISOitONIlItNCTIS) 

SUNl>SUN2sB*U 

NC1«NC2«0 

DO  ir  RIOz1*UMULT 

IFtRIRlOl  .LT<  RAV62>  60  TO  19 

SUM1<SUM1«A(K10>/S(K10I 

NC1«NC1«1 

60  TO  2' 

IS  CONTINUE 

SUH2sSUN2«ACKia>/SIKlC) 

NC2>NC2*1 
20  CONTINUE 
IS  CONTINUE 
NCTlIIzNCl 
NCT(2lzNC2 
OMIlIzSUHl/NCl 
0NI2IZSUN2/NC2 
RETURN 
END 
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I 

1 


I 


1 


5 


10 


15 


2U 


25 


30 


35 


40 


46 


5 . 


SUBOOUTINC  CNCR5T<NmiLT«R«V52«A«B«llteNtNCT» 

C SUBROUTIHC  CNCRBT  CONTUTCS  THC  RVCRRBC  Of  THC  OfSTRNCC  BCRSURCS 

C LNlHt  LRIL*  LUH(  AND  LH2L  FOR  THC  HORO  BCIMB  FROCCSUO. 

OIHCNSION  AI3»*2>«BI3e*2>*RC3e>*0NI4>*NCTI2l 

RL1R1SB.B2 

RLIH2«2.4A 

sun 1 >sun2  «sun3«sun4  • t 

NClzNC2«e 

00  1C  R10«1«NRULT 

C TIC  RCSIOUAL  DISTANCC  RCASURCS  ARC  LINCARIZCO  FOR  CACH  FRAHC. 
no  5 K5zl«2 

IFIB<KlStK5>  .LC*  RLinZI  60  TO  25 
RH2zB(Nl««R5l 
RHls(Rn2-RLin2*1.0l»RLini 
60  TO  30 
25  RnisRLini 
Rn2zRLtn2 

30  AA>1.0/<Rni-Rn2l 
BB>-AA*Rn2 

]FlAINlCtK5>  «LT.  RRII  60  TO  35 
IF(AIKl;*K5>  .6T.  RH2I  60  TO  4ft 
AIK10tK5lsAA*AIK10*K5t*BB 
60  TO  45 

33  A(N13«R5|sl.S 
60  TO  45 

40  t(R10«R5>s|..0 
45  CONTINUe 
3 CONTINUC 

C THC  H16H  AttO  LOtf  CHCRST  FRAHCS  ARC  FOUNO  ANO  THC  LINCAR  OISTANCC 

C nCASURFS  LHIH*  LHIL*  LH2H*  ANO  LH2L  ARC  COHAUTCO. 

IF«RIK1S>  .LT.  RAV62>  60  TO  15 

SUHlsSUni«A(NlC«ll 

SUn2sSUn2«A«K10«2> 

NClzNCl*! 

60  TO  ij 
15  CONTINUC 

sun3xsun3«A iNirtii 
SUn4:SUn4*A<Nlfat2> 

NC2:NC2«1 
20  CONTINUC 

1'  continue 

C the  4 AVCRAGC  OISTANCC  nCASURCS  FOP  THE  NOAO  BC1N6  ABOCESSCO  ARC  COHPUTCO. 
NCTI IISNCI 
NCT(2>sNC2 

IFCNCl  .C3.  41  NClsl 

1FINC2  .CO*  I NCZsl 

ON(l>s$!iMl/NCl 

on<2>zsun2/Nci 

On« J»sS0n3/NC2 

0NI4>z$un4/NCk 

PfTURN 

fcNC 


i 

1 


1 

i 


h 
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